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Abstract 

Interactions between the octet-baryons {B^) in the spin- flavor SUq quark model are inves- 
[ tigated in a unified coupled-channels framework of the resonating-group method (RGM). The 

interaction Hamiltonian for quarks consists of the phenomenological confinement potential, the 
color Fermi-Breit interaction with explicit flavor-symmetry breaking (FSB), and effective- meson 
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I exchange potentials of scalar-, pseudoscalar- and vector-meson types. The model parameters are 

' determined to reproduce the properties of the nucleon-nucleon (NN) system and the low-energy 

cross section data for the hyperon-nucleon interactions. Mainly due to the introduction of the 
vector mesons, the A^A^ phase shifts at non-relativistic energies up to Tiab = 350 MeV are greatly 
improved in comparison with the previous quark-model A^A^ interactions. The deuteron properties 
^ ' and the low-energy observables of the BgBg interactions, including the inelastic capture ratio at 

P5 ! rest for the scattering, are examined in the particle basis with the pion-Coulomb correction. 

The nuclear saturation properties and the single-particle (s.p.) potentials of Bg in nuclear medium 
are examined through the G-matrix calculations, using the quark-exchange kernel. The S s.p. po- 
tential is weakly repulsive in symmetric nuclear matter. The s.p. spin-orbit strength for A is very 
I small, due to the strong antisymmetric spin-orbit force generated from the Fermi-Breit interaction. 

The qualitative behavior of the BgBg interactions is systematically understood by 1) the spin-flavor 
SUq symmetry of Bg, 2) the special role of the pion exchange, and 3) the FSB of the underlying 
quark Hamiltonian. In particular, the B^Bg interaction becomes less attractive according to the 
increase of strangeness, implying that there exists no B^Bg di-baryon bound state except for the 
deuteron. The strong AA^-SA^ coupling results from the important tensor component of the one- 
pion exchange. The AA-SA'^-SS coupling in the strangeness S = —2 and isospin / = channel is 
relatively weak, since this coupling is caused by the strangeness exchange. The BgBg interactions 
are then applied to some of few-baryon systems and light A-hypernuclei in a three-cluster Faddeev 
formalism using two-cluster RGM kernels. An application to the three- nucleon system shows that 
the quark-model A^A^ interaction can give a sufficient triton binding energy with little room for 
the three-nucleon force. The hypertriton Faddeev calculation indicates that the attraction of the 
AA^ interaction in the ^Sq state is only slightly more attractive than that in the state. In the 
application to the aaA system, the energy spectrum of ^Be is well reproduced using the aa RGM 
kernel. The very small spin-orbit splitting of the ^Be excited states is also discussed. In the AAa 
Faddeev calculation, the NAGARA event for ^y^He is found to be consistent with the quark-model 
A A interaction. 
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1 Introduction 

Since the discovery of quantum chromodynamics (QCD), many efforts have been made to understand 
the structure of hadrons and hadron-hadron interactions based on this fundamental theory of the strong 
interaction. See, e.g., a review article [1] for these efforts. Such studies are of crucial importance in 
understanding many-body systems of nuclei and hadronic matter. A direct derivation of these elements 
from QCD, however, involves tremendous difficulties such as the quark confinement and multi-gluon 
effects in the low-energy phenomena. This has motivated a study using other effective models which take 
account of some important QCD characteristics like the color degree of freedom, one-gluon exchange, 
etc., even though they may not be rigorously derived from QCD. 

In this review article, we will focus on the interactions between the octet-baryons (Sg), studied in 
the non-relativistic quark model (QM), and their applications to light nuclear systems. The QM study 
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of the baryon-baryon interaction started from applying the three-quark (3g) model for the nucleon (A^) 
to the (3g)-(3g) system, in order to obtain a microscopic understanding of the short-range repulsion of 
the nucleon-nucleon (NN) interaction. Almost 30 years ago Liberman [2] first calculated the quark- 
exchange kernel of the phenomenological confinement potential and the (fTj • <Ty)(Af • A^)-type color- 
magnetic interaction originating from the Fermi-Breit (FB) interaction. The latter was suggested by 
De Rujula, Georgi and Glashow [3] as a possible candidate of the quark-quark (qq) residual interaction 
based on the asymptotic freedom of QCD. The effect of antisymmetrization among quarks is most 
straightforwardly taken into account in the framework of a microscopic nuclear cluster model, the 
resonating-group method (RGM). The RGM equation for the {3q)-{3q) system was solved by Oka and 
Yazaki [4] not only for the A^A^ interaction, but also for the AA interaction. Harvey [5] pointed out 
the appearance of the so-called hidden-color components in the many-quark systems. The nature of 
the repulsive core of the NN interaction was studied by Faessler, Fernandez, Liibeck and Shimizu [6] 
in the framework of RGM, including NN, AA and the hidden-color channels. These early studies 
were followed by a number of works [7, 8, 9, 10, 11], which repeatedly confirmed that the medium- 
and long-range effects, usually described by meson-exchange potentials, are absent in the QM and 
must be incorporated to obtain a realistic description of the baryon-baryon interaction. The first such 
attempt for the NN interaction was carried out [12] by supplementing the well-established one-pion 
exchange potential (OPEP) and a phenomenological medium-range central potential in the Schrodinger- 
type equation equivalent to the RGM equation. The deuteron properties were investigated in relation 
with the compact (6q')-bag configuration [13]. The most successful calculation of the NN interaction 
in this approach was given in Rcf. [14]. The authors introduced the spin-spin and non-central terms 
of the OPEP, assuming that the pions directly couple with the quarks in the (3g) core, and thereby 
yielding an effective-meson exchange potential (EMEP) at the baryon level. An alternative approach 
is to assume OPEP between the quarks in the RGM formalism and to calculate the quark-exchange 
kernel explicitly [15]. This program was carried out by the Tubingen group [16, 17], who extended the 
model of the A^A^ interaction to the hyperon-nucleon {YN) and hyperon-hyperon {YY) interactions as 
well. Their approach includes a pseudoscalar (PS) meson exchange between the quarks as well as a 
phenomenological cr-like potential at the baryon level. A fully microscopic calculation including both 
PS- and cr-meson exchange potentials at the quark level was undertaken by the Salamanca group [18, 19] 
for the A'^A'^ interaction in the framework of the chiral QM. This model was further extended by the 
Beijing group [20, 21, 22, 23] to a chiral SU^ model with strangeness degree of freedom, and applied to 
the YN and YY interactions. In the meantime, meson-exchange models have also made a progress in 
the description of the medium-range attraction of the YN interaction. The new Jiilich YN model [24] 
has discussed the contributions of the scalar-isoscalar (a) and vector-isovector (p) exchange channels, 
constrained by a microscopic model of correlated nn and KK exchanges. The extended Nijmegen 
soft-core model ESC04 [25, 26] has introduced two PS-meson exchange and meson-pair exchange, in 
addition to the standard one-boson exchange and short-range diffractive exchange potentials. 

Our viewpoint for applying the RGM formalism to the baryon-baryon interactions is that the SUq 
QM with some realistic modifications like peripheral mesonic or (qq) effects can essentially reproduce the 
whole low-energy hadronic phenomena. For the low-energy phenomena related to the S-wave ground- 
state hadrons, the predictions with the FB interaction of the one-gluon exchange have no apparent 
disagreement with experiment except for some trivial cases. Though the FB spin-orbit (LS) interaction 
is believed not to be favorable for the very small energy splitting of the negative-parity nucleon excited 
states, known as the "missing LS force problem" of the P-wave baryons, it was argued in Ref. [27] 
that this problem can be resolved by a proper treatment of a compact (3g) structure embedded in 
the (3g)(gg) continua. In other papers [28, 29], we showed that the same FB LS force has almost 
right order of magnitude compatible with the experimental information on the A^A^ and YN spin- 
orbit interaction. Fujiwara and Hecht discussed in [30, 31] the off-shell FB interaction in the non- 
relativistic framework and incorporated the meson-exchange effects to the {3q)-{3q) RGM by exphcitly 
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introducing (gg) and (gg)^ excitations. As the main feature of the A^A^ interaction was reproduced in 
the extended QM without introducing any extra meson parameters, it may be considered the first step 
toward the conventional meson-exchange description through a possible mechanism of {qq) exchanges 
between nucleons. An extension of the model to include the antiquark degree of freedom is also made 
to describe other hadronic interactions such as the pion-pion and pion-nucleon scatterings as well as 
the proton-antiproton annihilation processes [32, 33, 34]. 

We note that to apply the RGM formalism to the (3g)-(3g) system is not straightforward, since there 
are some new features arising from the intrinsically relativistic nature of the system. The standard non- 
relativistic RGM formulation should be augmented with the following features: 1) the full FB interaction 
with flavor-symmetry breaking (FSB), including the momentum-dependent Breit retardation term and 
the Galilean non-invariant LS term, 2) a (Os)^ harmonic-oscillator (h.o.) wave function with a common 
size parameter b for all the (3g) systems, 3) a perturbative treatment of the FSB, in terms of the pure 
SUq QM wave functions with the [3]-symmetric spatial function, and 4) correct reduced masses and 
threshold energies in the RGM formahsm. See Refs. [35, 36] for the background of these extensions. 

The effective qq interaction in our model consists of a phenomenological quark-confining poten- 
tial, the one-gluon exchange FB interaction, and other terms generating EMEP's from various meson- 
exchange mechanisms. Various versions of our QM for the baryon-baryon interactions emerge from 
several possibilities of how to incorporate these EMEP's in the RGM equation. A simplest method is to 
use effective local potentials V^^ at the baryon level, which are determined from the existing one-boson 
exchange potentials (OBEP) with some adjustable parameters. From the correspondence between the 
RGM formalism for composite particles and the Schrodinger equation for structureless particles, the 
relative-motion wave function should be rcnormalized with the square root of the normalization kernel 
\/N. a nice point of this prescription is that we only need the normalization kernel to generate the ex- 
change terms. Our earhest version of the QM, RGM-F [35, 36, 37], uses this prescription for EMEP's, 
where V^^ is assumed to be proportional to the scalar- (S-) meson and PS-meson exchange potentials 
of the Nijmegen model-F [38]. More specifically, the central force of the S-meson nonet and the tensor 
force of the tt and K mesons in the Nijmegen model-F are used. With only two adjustable parameters 
determined in the A^A^ sector, RGM-F can reproduce reasonably well all the low-energy cross section 
data of the YN systems. 

In the next step, we have improved the way of introducing the EMEP's in two respects. One is 
to calculate the spin-flavor factors exactly at the quark level and the other is to include the spin-spin 
terms originating from all the PS mesons. We have shown in Ref. [39] that the available A^A^ and YN 
data can simultaneously be reproduced in the standard (3g)-(3g) formulation, if one assumes the full 
PS- and S-meson nonet exchanges at the quark level and properly introduces the FSB in the quark 
sector. An explicit evaluation of quark-exchange kernels for the EMEP leads to a strong constraint 
on the flavor dependence of each A^A^ and YN channel. The SU^ relation of the coupling constants 
emerges as a natural consequence of the SUq QM. For S mesons, one of the SU3 parameters, F/{F+D) 
ratio, turns out to take the SUa value of a purely electric type. Under this restriction, the low-energy 
A" A", E+p and Ap total cross sections cannot be reproduced consistently with a unique mixing angle of 
the flavor-singlet and octet S mesons. We resolve this problem in two ways; one is to change the mixing 
angle only for S+p channel and the other is to employ the same approximation as RGM-F, changing 
the F/{F + D) ratio, solely for the isoscalar S mesons, e and S*. These models are called FSS and 
RGM-H [39, 40, 41], respectively. Since their predictions are not much different except for the roles 
of the LS^^'' force in the AN-T,N{I = 1/2) coupled-channels (CC) system, the result with FSS will be 
discussed in this article. 

In the most recent model called fss2 [42, 43], we have also included the following interaction pieces 
mediated by S and vector (V) mesons: LS, quadratic spin-orbit (QLS), and momentum-dependent 
Bryan-Scott terms. The introduction of these terms to the EMEP's is primarily motivated by the 
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insufficient description of the NN experimental data with our previous models. Due to these improve- 
ments, the A^A^ phase shifts at non-relativistic energies up to Tiab = 350 MeV have attained an accuracy 
almost comparable to that of OBEP's. The model parameters are searched in the isospin basis to fit 
the A^A^ phase-shift analysis with J < 2 and Tiab < 350 MeV, the deuteron binding energy, the ^So 
NN scattering length, and the low-energy YN total cross section data. The model fss2 reproduces the 
YN scattering data and the essential features of the AA^-SA^ coupling equally well with our previous 
models. Another important progress for the QM baryon-baryon interaction has come from a technique 
to solve the RGM equation in the momentum representation. Here we solve the T-matrix equation 
of the RGM kernel, for which all the Born terms are derived analytically. The partial wave decompo- 
sition of the Born kernel is carried out numerically in the Gauss-Legendre integral quadrature. This 
method, called Lippmann-Schwinger RGM (LS-RGM) [44], enables us not only to obtain very stable 
and accurate results for the ^-matrix over a wide energy range up to a few GeV, but also to apply our 
QM baryon-baryon interactions to physically interesting problems such as G-matrix and three-cluster 
Faddeev calculations. The single-particle (s.p.) potentials of the Bg (A", A and E) and the strength of 
the s.p. spin-orbit potential are predicted through the G-matrix calculations [45, 46]. 

We have carried out coupled-channels RGM (CCRGM) calculations both in the isospin basis and in 
the physical particle basis. The former is convenient to show concisely the characteristics of the BsBg 
interactions, but the latter is necessary to make a realistic evaluation of cross sections for a comparison 
with experimental data. The latter is also necessary for the study of the charge dependence of the A" A" 
interaction and the charge-symmetry breaking (CSB) of the AA^ interaction. The charge-independence 
breaking (GIB) is partially explained by the so-called pion-Goulomb correction [47], which implies 1) 
the small mass difi^crence of the neutron and the proton, 2) the mass difference of the charged pions 
and the neutral pion, and 3) the Goulomb effect. For the GSB of the AA^ interaction, the early version 
of the Nijmegen potential [48] already focused on this GSB in the OBEP including the pion-Goulomb 
correction and the correct threshold energies of the AN-YIN coupling in the particle basis. We always 
include the pion-Goulomb correction when the particle basis is used for the calculation. For the Goulomb 
effect, we calculate the full exchange kernel without any approximation. The pion-Goulomb correction 
and the correct treatment of the threshold energies are found to be very important to analyze the 
low-energy observables in the HN-AN GG problem. 

In the next section, we briefiy describe the formulation of the (3g)-(3g) RGM and its application to 
the G-matrix and Faddeev calculations. These formulations are made in the momentum representation. 
The essential points to be emphasized include the energy dependence of the RGM kernel and the exis- 
tence of the Pauli-forbidden state. The treatment of the Goulomb force in the momentum representation 
is also given. Section 3 presents results and discussions. We start with the two-nucleon properties and 
proceed to the YN interaction, the G-matrix calculations, the YY interaction and the other BgBg in- 
teractions including the SN, SY and SS interactions. The importance of the spin-flavor SUq symmetry 
in our QM framework is stressed through a systematic analysis of the phase-shift behavior and total 
cross sections in some particular channels. In Sec. 3.4, our QM baryon-baryon interactions, FSS and 
fss2, are applied to the calculations of the triton and hypertriton binding energies in the three-cluster 
Faddeev formalism. The application to other light hypernuclei, ^Be and aaHc, is given in Sec. 3.5. The 
spin-orbit splitting of ^Be excited states is also examined by the Faddeev and G-matrix calculations. 
A summary is drawn in Sec. 4. 
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2 Formulation 



2.1 RGM formalism for quark clusters 

The RGM is a microscopic theory which describes the dynamical motion of clusters, on the basis of 
two- (and three-, • • • ) body interactions of constituent particles. For the application to baryon-baryon 
interactions in the QM, the Sg's are assumed to be simple s-shell clusters composed of the (Os)^ h.o. 
wave functions with a common size parameter b. They have spin-flavor and color degrees of freedom as 
well, for which we follow the standard non-relativistic SUq description. We use the Elliott notation to 
specify the SU3 quantum numbers of Bs, and write Bg = -B(ii)a with a = YI (hypercharge Y, isospin 
/) to denote a flavor function N {YI = 1 i), A (YI = 0), S (YI = 1) or S (YI = -1 1). We also 
use the notation a — alz to denote a specific member of the isospin multiplet for each baryon. Then 
the QM description of the Bs is achieved with the correspondence 

B^n)a - W^I^^JUS) = YI 4f [[^i(lVi(2)]5m(3)]i 

^ S'(A^)=0(01), 1(20) ^ 

X [ [ F(io)(l)F(io)(2) F(io)(3) ](ii)a , (2.1) 

because the spin-flavor part of B^ is totally symmetric. Here wi is the spin wave function and -F(io) the 
flavor wave function of a single quark. The QM wave functions of the Bs are expressed as 

0^3(123) = 0(-b)(123) W[%^,,^^{123) C(123) , (2.2) 

where C(123) is the color-singlet wave function. The orbital wave function, (f)^°^^\l23) , is an intrinsic 
part of Bg and contains no dependence on the cm. coordinate Xq — {xi + X2 + Xs)/^. The orbital 
functions for the (3Q')-clusters are thus assumed to be flavor-independent and are taken to be the same 
for aU the -Bs's. 

The RGM wave function for the (3g)-(3g) system can be expressed as 

* = J]^'{te(i?)} , (2.3) 

a 

where the channel wave function 0q, = 0*^°"^^^^^^^*^ is composed of the orbital part of the internal 
wave function 0^°'''^) = 0^°'''^) (123)0*^°'''^-' (456), the isospin-coupled basis of the spin-flavor SUq wave 
functions and C,^ = C(123)C(456). A detailed definition of C,a^ with a set of specific two-baryon 
quantum numbers a will be given below. The antisymmetrization operator A' in Eq. (2.3) makes ^ 
totally antisymmetric under the exchange of any quark pairs and can be reduced to the form A' — > 
(1/2)(1 — 9P36)(1 — Pq) with the core-exchange operator of the two (3g)-clusters Pq — -Pi4-f25-P36- 

We now explain the basis . In the general B^Bg interaction, there appear several new features 
which are absent in the A^A^ interaction. These are mainly related to the rich fiavor contents of the BgBg 
system constrained by the spin-flavor SUe symmetry and the generalized Pauli principle. Since these 
features are common in both OBEP and QM descriptions, we will use the baryon's degree of freedom 
as much as possible. Here we use the isospin-coupled basis for presentation, but the transformation to 
the physical particle basis is straightforward. The -SC/s-coupled basis of two Ps's is defined by 

[B{n)B(u)] = [5(11) (1) 5(11) (2)] 

= J]((ll)ai(ll)a2 II (A/x)a)45 (ll)ai(l) -B(ii)a2(2)]^^^ , (2.4) 

0-1,0,2 

where the last brackets denote the isospin coupling and p denotes a possible multiplicity label of the 
SU3 coupling, (11) X (11) (22) + (30) + (03) + (11), + (11)„ + (00) (p = 1 for (11)„ and p = 2 for 
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(ll)s). We assume that the first i?(ii) is always for particle 1 and the second for particle 2. A set of 
the internal quantum numbers, a {a = alz), stand for the hypercharge Y = Yi + Y2, the total isospin / 
and its z-component J^, which are all conserved quantities if we assume the isospin symmetry and the 
charge conservation Q/e — Y/2 + 1^. Only when the particle-basis is used, the isospin / is no longer 
"exactly" conserved. The flavor-symmetry phase V of the two-baryon system is determined from the 
symmetry of the SU^ Clebsch-Gordan coefficients through 

P12 [^(11)^(11)] (A^)a;p = ^ [^(11)^(11)] (A^)a;p ' (2-5) 



with 



Here P^^ 



(— 1)"^+^ in multiplicity- free case , 
(-1)^ for (A/x) = (11) . 

(Af ■ Af)/2 + 1/3 is the ffavor- exchange operator with Af 



(2.6) 

1, 2) being the Gell-Mann 



matrices for the ffavor degree of freedom. It is important to maintain the deffnitc ffavor symmetry even 
in the isospin-coupled basis and in the particle basis, since the ffavor symmetry is related to the total 
spin and the parity through the generalized Pauli principle; 

\L+S'p 



-1 



1 



(2.7) 



where L is the relative orbital angular momentum. For the NN system with Y = 2, the total isospin 
/ is good enough to specify V; i.e., V = (-1)^"^ so that Eq. (2.7) becomes the well-known rule 
(—1)-'""'"'^"''^ — —1. More generally, the flavor-exchange symmetry V is actually redundant for the systems 
of identical particles with ai = 02 and Ji = I2, since it is uniquely determined as P = (— 1)^^^"''^. Let 
us introduce a simplified notation Bi = -B(ii)ai {i — l,---,4 in general) with Oj = y^/j. Then the 
symmetrized isospin basis is defined through 



[BiBil^j 




J2 ((ll)ai(ll)a2 II (A/i)a;p) 



(2.8) 



(A/i)a;p 

where the sum is over {XfJ.)p compatible with Eq. (2.6). As is apparent from the construction, the basis 

in Eq. (2.8) is naturally an eigenstate of the fiavor-exchange operator with the eigenvalue V. It is 
now straightforward to incorporate the spin degree of freedom in the above formulation and extend the 
argument to the spin-fiavor SUq wave functions of the {3q)-{3q) system. In the QM, the isospin basis 
with a definite flavor-exchange symmetry is given by 

1 



SF 



(ll)ai 



[3] 

2(11)02^ 



+ n-1) 



h+h-I 



^i(n)a.(123)VF|;;^, (456) 



(ll)a2 



[3] 

2(ll)ai' 



The subscript a in Eq. (2.9) specifles a set of quantum numbers of the channel wave function; 

SS,YIh-V . 



a 



^(11)«4(11)«2 



(2.9) 



(2.10) 



In the particle basis, we use the notation Bo^ 



Oi2J 



^J2{1 (^ai.aa) 



B{ii)ai s-^d deflne 

{BaiBa2 + VBa-^Bai] 



■'a\,a2 



1 + 5, 



01,02 J J 



V 



(2.11) 
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The channel quantum number a in Eq. (2.10) in this case is specified by a = [|(11) ai|(ll)a2] SSz', V 
with Y = Yi + Y2 and = Ji^ + l2z- For the system of identical particles with cti = 0:2, the flavor- 
exchange symmetry is again redundant and V — 1. The relationship between the isospin basis and the 
particle basis in Eq. (2.11) is inversely written as 



hzhz 



01,02 



The bases in Eqs. (2.9) and (2.11) arc the cigenstates of the core exchange operator Pq = Pf/ P25 -^36 
with the eigenvalue {—lY~^V, where P,^j^ = PijPfj exchanges the spin-flavor coordinates. 

The QM Hamiltonian consists of the rest-mass plus non-relativistic kinetic-energy term, the quadratic 
confinement potential, the full FB interaction with explicit quark-mass dependence, and the S-, PS-, 
and V-meson exchange potentials acting between quarks: 



H 



t ("•.-' + ir.-Ta)+t + Uf- + j: Uff + E C/-^ + E U:/] . (2.13) 



Here U^^ = — (Af ■ A^^) Qct'^ with r = \r\ = \xi — Xj \ is the confinement potential of quadratic power law, 
which gives a vanishing contribution to the interaction in the present formalism. The qq FB interaction 
C/™ is composed of the following pieces [7, 49] : 

-Ug^ + Ulf^ + Ug^ + U:l^' + U^' + , (2.14) 

where the superscript (CC) stands for the color-Coulombic or (Af ■A^)/r piece, (MC) for the momentum- 
dependent Brcit retardation term or (Ap • A^)|(j>j ■ Pj) + -Pi) 'Pj/f^^^ / {mimjr) piece, {GC) for the 

combined color-delta and color-magnetic or (Ap ■ A^) |l/(2m^) -|- l/(2mp -|- 2/(3mjmj)(crj ■ crj)|(5(r) 

piece, {sLS) for the symmetric LS, (aLS) for the antisymmetric LS, and (T) for the tensor term. It 
is important to note that all the contributions from the FB interaction are generated from the quark- 
exchange diagrams, since we assume color-singlet cluster wave functions. The EMEP's, [/^^, t/^."'^, and 

C/j^^ will be discussed in the next section. When the calculations are made in the particle basis, the 
Coulomb force is also introduced at the quark level. 

The RGM equation for the parity-projected relative wave function Xq(-R) is derived from the vari- 
ational principle — = 0, and it reads as [11, 40] 

'"'^^a (m) ^^^^^ " ? / R', E) xl'W , (2.15) 

where M. aa' {R, B! ; E) is composed of various pieces of the interaction kernels as well as the direct 
potentials of EMEP: 

Mo^AR. R, E) = 5{R -R')J2Y1 vSL'iR) + E ^aa'iR, R') - M^AR, R) ■ (2-i6) 

The relative energy in channel a is related to the total energy E of the system in the cm. system 
through Ea = E — E"^^, where i?™* = E"^^ + i?™* with a = aia2. In Eq. (2.16) the sum over Q in the 
direct term implies various contributions of interaction types for the meson-exchange potentials, while 
(3 specifies the meson species. These EMEP interaction types are classified to i7 = CN, SS, TN, LS, 
and QLS, which are the spin-independent central, spin-spin, tensor, LS, and quadratic LS terms. 
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respectively. On the other hand, Q for the exchange kernel Ai^^,{R, R') involves not only the exchange 
kinetic-energy (K) term but also various pieces (CC, MC, GC, sLS, aLS, T) of the FB interaction, 
as well as the EMEP contributions corresponding to the direct term. The structure of the EMEP 
contributions is rather involved and contains some ambiguities in the way of introducing them even 
starting from the meson-exchange potentials at the quark level, as will be discussed in the next section. 
Here we only refer to a general procedure to solve the RGM equation (2.15) and a couple of important 
features of the present RGM framework using solely qq interactions. 

A simple way to solve the CCRGM equation (2.15) is the variational method developed by Kamimura 
[50]. Although this technique gives accurate results up to about Tiab = 300 MeV, it seems almost 
inaccessible to higher energies due to the rapid oscillation of the relative wave functions. We alternatively 
use a method to solve the CCRGM equation in the momentum representation, namely, an LS-RGM 
equation [44]. The idea of solving the RGM equation in the momentum representation is not new [51]. It 
is also used by Salamanca group in the study of the QM baryon-baryon interactions [52] . In this method 
all the necessary Born amplitudes (or the Born kernels) for the quark-exchange kernels are analytically 
derived by using a transformation formula, which is specifically developed for momentum-dependent 
two-body interactions acting between quarks [44]. The partial- wave decomposition of the Born kernel 
is carried out numerically in the Gauss- Legendre integral quadrature. The LS-RGM equation is then 
solved by using the techniques developed by Noyes [53] and Kowalski [54] . Although this method requires 
more CPU time than the variational method, it gives very stable and accurate results in a wide energy 
range. Since we first calculate the Born amplitudes of the RGM kernel, it is almost straightforward 
to proceed to the G-matrix calculations [45], as well as many Faddccv calculations in the three-cluster 
Faddeev formalism. The transformation to the RGM equation (2.15) is easily carried out by simple 
Fourier transformations. 

A nice property of the present RGM formalism is that the internal-energy contribution is already 
subtracted in the exchange kernel. Namely, A4^^^,[R, R') in Eq. (2.16) for the central components 
O, — K, CC, MC, GC, CN, and SS is defined through its corresponding original exchange kernel 
Mf^r^iR^R'): 

Mf^,(R, R) = Mf]^-\R, R!) - + <)) M^R, R') , (2-17) 

where denotes the i7-term contribution to the internal energy of the Bg specified by the fiavor 

label a = 0102. Owing to this subtraction, the mass term of the kinetic-energy operator in Eq. (2.13) 
as well as the confinement potential with Q = Cf exactly cancels between the first and the second 
terms on the right-hand side of Eq. (2.17). We consider this feature one of the advantages of the RGM 
formalism, because the present QM produces those results which are independent of the strength of the 
confinement potential and are insensitive to the details of the confinement phenomenology. 

We should note that the original Hamiltonian (2.13) in general contains the derivative operators 
d/dxi and d/dxj, which cannot be reduced to the derivative of the relative coordinate r — Xi — Xj be- 
tween two quarks. The appearance of these Galilean non-invariant terms, like the momentum-dependent 
Brcit retardation term U^'" and the antisymmetric LS term U°j"^ , is a direct consequence of the more 
strict Lorcntz invariance at the relativistic level, and their RGM kernel should be explicitly evaluated 
in the total cm. system [28, 35]. We include the U^'-^ term and take account of its contribution to the 
relative kinetic-energy term explicitly. Even with this procedure, the calculated reduced mass of the TiN 
system is degenerate with that of the AA^ system. This is one of the limitations of the non-relativistic 
QM, in which the inertia masses (of E and A in the present case) are not always reproduced correctly. 
In order to use the correct reduced masses in the CCRGM equation, we make the following replacement 
only for the exchange kernels M.l^}{R, R') and Ail^F\R, R'): 



Mf^,{R,R!)^Mf^,{R,R!)^^,Mf^,{R,R) Q^K and MC , (2.18) 



exp 
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where «o is the incident channel. (This prescription is shghtly different from that of Refs. [37, 40].) For 
such a case that the exact Pauh-forbidden state exists in the CC system, a modification is necessary as 
discussed in Sec. 4 of Ref. [55]. The essential idea is as follows. Let us express the CCRGM equation 
(2.15) as 

(e-Ho- l^^°^)x = , (2.19) 

with e = E-Ei^t and V^^^ = Vi) + G + eK. Here V^^^ represents the RGM kernel in Eq. (2.16) with 
Vd the direct potential, G the exchange kernel '^^Ai^^,{R,R'), and K the exchange normalization 
kernel —Ai^^,{R, R'). The essential feature of the RGM equation (2.19) is the existence of the trivial 
solution X = u, which satisfies Ku = u. If we use the projection operator on the Pauli-allowed space, 
A = 1 — \u){u\, Eq. (2.19) can equivalently be written as 

A(£ -Hq- V^^^)Kx = . (2.20) 

In this form, we can safely replace the calculated reduced mass with the empirical mass ix^^"^ . Quite 
obviously, this procedure is accurate only when the calculated mass is a good approximation for the 
empirical mass. The same procedure is also applied when we formulate the CCRGM equation in the 
particle basis. The present constituent QM does not reproduce the small mass difference of the isospin 
multiplets exactly even if the FSB of the quark masses and the Coulomb interaction at the quark level 
are properly introduced. We therefore solve the modified equation 

^^^exp _ ^exp _ ^RGM-)^^ ^ q ^ ^2.21) 

instead of Eq. (2.20), where £^^p and //q^'p are defined by 

£-P ^E- E^l , Hi'''' = — ( = -^Ho . (2.22) 

This procedure to use the empirical reduced masses and threshold energies in the Pauli-allowed model 
space has a wide applicability to remedy the inflexibility of the RGM framework. If we restore Eq. 
(2.22) to the original Schodinger-type equation (2.19), we obtain an additional correction term to the 
RGM kernel, which is supposed to be small enough for practical calculations, but necessary to preserve 
the exact treatment of the Pauli principle: 

(s^'^P - HI'''' - V^^^ - AG)x = . (2.23) 

Here AG is given by 

AG = A {AEi^, + AHo) A - {AE,^, + AHo) (2.24) 

with 

A^;^., = E^:^ - E,^, , AHo = -Ho^ (^-^ - 1^ Ho . (2.25) 

A detailed example is given in Ref. [55] for the AN-T,N[I = 1/2) CCRGM, which involves a Pauli- 
forbidden state in the (ll)s SU^ representation for the 'So state [37]. With this procedure we can 
employ the empirical reduced masses and threshold energies without impairing the major role of the 
Pauli principle. 
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2.2 Effect ive-meson exchange potentials 



In this subsection, we outline a procedure to derive the RGM kernel, starting from a general form of 
the qq interaction, and discuss appropriate forms of the EMEP's used in the model FSS and fss2. Refer 
to the original papers [56, 57, 40, 42, 44] for detailed calculations. 

The systematic evaluation of the quark-exchange kernel is carried out for a two-body interaction 

C/g = a"<^g < , (2.26) 

where is a simple numerical factor including a single quark-gluon coupling constants as — g^d (47r) 
and various types of quark-meson coupling constants Qqq^-, fqq^, etc., w^j represents the spin-flavor part, 
wfj the color part and ufj the spatial part. The color part is wfj = (l/4)(Af ■ A^) for the FB interaction 
and wfj = 1 for EMEP's. In Table 1, a^, wfj and ufj are given for the FB interaction [7], while 
those for EMEP's in Tables VII and VIII of Ref. [42]. Four different types of the spin-fiavor factors 
fl — C, SS, T, LS are required for the most general EMEP up to the V mesons; wfj = 1, w^j^ = {ai-crj), 

"^Ij — '^if — (^i + For the flavor-octet mesons, these spin operators should be 

multiphed with (Af • Af ). The non-central factors are defined by the reduced matrix elements for the 
tensor operators of rank 1 and 2. For example, the tensor operator is expressed as 

(fe, k) = 3 {(Ti ■ k) {(Tj ■ k) - {cTi ■ (Tj) k^ = sVlQ 

with y2^J,{k) = ^^/AttJTE k^Y2^{k) . The Born kernel for the two-body interaction (2.26) is calculated by 

6 

M{qf, q,) = { e'^r^ | Ur,A^ I 4>) = Y. M^q^, q,) . (2.28) 

i<j xT 

where (f) = (f)^°^'°^ ^ with ^ = is the h.o. internal wave function of the (3g)-(3g) system and the 

superscript ^2 for specifying the interaction piece, as well as the channel subscript a, are omitted for 
simplicity. In Eq. (2.28), the spatial part of the Born kernel is expressed as 

M,r(q/,q,) = (^.e^^/-^0(°'^^)|ii,,|e^«^-^0(°'^'^)) with (i,j)eT. (2.29) 



1(2) 



y2{k) 



(2.27) 



Table 1: The coefficients a^, the corresponding spatial functions u^{k,p^,p2), and the spin-flavor 
operators for the Fermi-Breit interaction. The quark pair i, j = 1, 2 is assumed. For the non- 
central terms, the angular-momentum coupling JJ-^ — q;^3\/I0 [m^w^] etc. should be taken. 



~~CC 

MC 

GC 



a^/as 



fc2 



1 

mud 



\my,d J 

sLS ff^V 
aLS -li^ 
T ^(-^\ 

12 \^ m„ri J 



[2, ■ p,) {k ■ p,) ~ {p, ■ p,)] 
1 

[k,l{p,-p,)] 

%i [fc, Pi + P2] 



1 

mim2 
2" 



m^d 
m2 



1 / rriud 
3 \ mi 



f 1 

2 



1 / rriud 
3 \ m2 



+ 



cr2 



2 (mud) 

3 mim2 



f "W^ (J -\- ( '"ltd 



cr2 + 



2 (rn^d) 

3 mim2 

m\m2 



(CTI • 0-2) 

(CTI + (T2) 
(cTi - 0-2) 



m\m2 



(2) 
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It is specified by the number of exclianged quarks, x = 0, 1, and various interaction types, T. More 
specifically, x = with Zq = 1 corresponds to the direct term (EMEP's only) and x = 1 with zi = P^q 
the (one-quark) exchange term. The interaction types, T = E, S, S', and D_, correspond to some 
specific (i, j) pairs of quarks [56]. The spin-flavor-color factor X^t is defined through 

r 

i<j 

where Co = 1 and Ci = —9. The color-factors X^j- = {zx \ w^j \ ^*") for each G T are given 
as follows. For the quark sector with w^j = (l/4)(Af ■ A^), = -(2/3), X,^^ = 0, and Xf^ = 
Xg = X^^s' = -(2/9), X^jy^ = 1/9, X^jj_ = 4/9. For the EMEP sector with wg = 1, X^^ = 1 and 
X^^ = 1/3. The spin-flavor-color factor for the exchange normalization kernel is defined by Xn — 
(-3)(P36 \^^^)- We also need Xk = 24(P36 | ^6 - ^5 | C^^) with the hypercharge operator Yi 
for the exchange kinetic-energy kernel of the YN systems. 

To calculate the matrix elements in Eq. (2.30), it is convenient first to calculate three-quark matrix 
elements of various operators with respect to the [3]-symmetric and [21] -symmetric SUq wave functions. 
We need a small number of operators for this calculation, which are always expressed as composite 
operators of the two SUq unit vectors [49, 58]. After representing the spin-flavor factors in {3q)-{3q) 
system as operator products of these composite operators, we can evaluate the (3g)-(3g) matrix elements 
of these operator products, using a simple formula which combines the (3g) matrix elements in various 
ways. The advantage of this method is that it is directly related to the straightforward separation of 
spin-isospin dependence of various non-central forces appearing in the BsBg interactions. The final 
results of some spin-fiavor factors are already published in our previous papers. In the quark sector, the 
operator form of the spin-fiavor factor is expressed as X^, which is essentially X^j- in Eq. (2.30) but 
has a slightly different numerical factor. The spin-fiavor-color factors in the quark sector for the NN 
system are given in Eq. (B.6) of Ref. [31]. We should note that Xj- there is actually X^'-'' for T = 
and and Xf^ = Xf,^ = -Xg. (Note that Xe is not Xf^.) Furthermore, X^^f = X^^ = Xn, 
X^^ = X|f^ = -2Xjv, X^^ = -Xn, and X^ = for A^A^. The spin-orbit and tensor factors for 
YN systems are found in Appendix C of Ref. [59]. The onc-quark exchange {x = 1) central factors 
{X^)B:iNBiN in the YN system with F = A, S and S are given in Appendix B.3 of Ref. [44]. 

We define the Born kernel for the exchange kernel Aiaa'{R, R ] E) in Eq. (2.16) as 

M^AQf. Qi: E) = ( e^«/-^ I Mo^AR, R, E) \ e^^> 

Each component of the Born kernel is given in terms of the transferred momentum k — and the 

local momentum q — {qf + qj)/2. In Eq. (2.31) the space-spin invariants — 0^{qf, q^) are given by 

Qcentral ^ ^ 

with n=[q,xq^], S'=i(o-i + <72) , 5(-) = i(o-i-o-2) , P, = i^t^^. (2.32) 
For the tensor and QLS parts, it is convenient to take four natural operators defined by 

- 5i2(/e, k) , O^' = S,2{q, q) , O^" = S.^ik, q) , = 5i2(n, n) , (2.33) 



for X 



(2.30) 
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where Si2{a,b) = (3/2) [ (cti ■ a) (erg ■ 6) + (0-2 ■ a) (cri ■ 6) ]- (cri ■ cts) (a ■ b). The calculation of Eq. (2.31) 
is most easily carried out by using a special type of the transformation formula, developed in Appendix 
A of Ref. [44] for the systems of two (Os) ^-clusters. This formula is particularly useful to calculate the 
Born kernel from the momentum-dependent two-body interaction. For the direct term of the static 
two-body interaction, this formula is reduced to the simple multiplication of the Gaussian form factor 
F{kY = e"'^^'^) to the two-body interaction in the momentum representation. For the exchange term, 
the exchange normalization kernel is factored out. The formula is especially simple for the Gaussian 
interaction, but is also applicable to the Yukawa functions preserving the_factorization property. The 
resultant Born kernel is expressed using generalized Dawson's integrals hn{x) and modified Yukawa 
functions ya{x), Za{x) etc., which arc given by the error function of the imaginary argument. The 
momentum dependence of the two-body interaction appears as the polynomial terms in q up to the 
second order. The direct Born kernel Af^|^^,(qj, q^) in the EMEP sector is shown below. The exchange 

Born kernel M^l{qj^, q^) is given in Appendix B of Ref. [44] for the FB interaction and in Appendix A 
of Ref. [42] for the EMEP's. 

The EMEP's at the quark level are most easily formulated in the momentum representation, by 
using the second-order perturbation theory with respect to the quark-baryon vertices. We employ 
the following qq interaction, which is obtained through the non-relativistic reduction of the one-boson 
exchange amplitudes in the parameter 7 = {m/2mud) (where m is the exchanged meson mass and mud 
is the up-down quark mass): 



An 



k +m'^ 



- 1 + 



in ■ S 



= -ff 



47r 



U 



Ait 



- (1 - c, 



k -\-m? 



k){<Tj • fe) - (1 - C5)(m^ + k^)-{(Ti ■ (Tj) 



1 + 



rf 



Imudm 



rriudm 



{(Ti ■ n){(Tj ■ n) 
-in ■ S 



(2.34) 



Here the quark- meson coupling constants are expressed in the operator form in the flavor space [58, 59]. 
For example, the product of the two different coupling-constant operators g and / are expressed as 



9f 



9ih 



for 



singlet mesons 
octet mesons 



9sfs^a^ai^aj 

The meson mixing between the flavor-singlet and octet mesons is very important, implying 

fr^> = ficos9 + fssinO Xs , fr, = -fisinO + fscos9 Xs , 



(2.35) 



(2.36) 



instead of /i and /sAg in Eq. (2.35) for the PS mesons. A similar transformation is also applied to 
the S- and V-meson coupling constants. The SU3 parameters of the EMEP coupling constants are 
therefore /i, /§ and 6. The S-meson exchange EMEP in Eq. (2.34) involves not only the attractive 
leading term, but also the momentum-dependent q^ term and the LS term. In the PS-meson exchange 
operator, the parameter cs is introduced only for the one-pion exchange, in order to reduce the very 
strong effect of the delta-function type contact term involved in the spin-spin interaction. The case 
Cs = 1 corresponds to the full expression, while Q = corresponds to the case with no spin-spin contact 
term. The V-meson exchange potential is composed of the electric-type term, the magnetic-type term 
and the cross term. In the electric term, the central force generated by the cu-meson exchange potential 
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is usually most important, and it also includes the g^-type momentum-dependent term. As to the 
introduction of the V-meson EMEP to the QM, the problem of double counting is addressed especially 
in the strong short-range repulsion from the uj meson [60]. We will not discuss this problem here, 
but avoid this double counting for the short-range repulsion and the LS force, by simply choosing 
appropriate couphng constants for V mesons, i.e., fY^, fY"^ ~ 1, and fg^ — (see Table 2). The 
magnetic term is usually important for the isovector p meson, and yields the spin-spin, tensor and QLS 
terms in the standard OBEP. The choice in Eq. (2.34) is to keep only the QLS term with the spin-spin 
term proportional to L"^, the reason for which is discussed below. Finally, the cross term between the 
electric and magnetic coupling constants leads to the LS force for the qq interaction. The antisymmetric 
LS (LS'(-)) force with S = (cri -cr2)/2 is not generated from EMEP, b ccause the flavor operator in Eq. 
(2.35) is the Gell-Mann matrix and also because the mass difference between the up-down and strange 
quarks is ignored in Eq. (2.34). 

We should keep in mind that these EMEP's are by no means the consequences of the real meson- 
exchange processes taking place between the quarks. First of all, the static approximation used to derive 
the meson-exchange potentials between the quarks is not permissible, since the masses of S mesons and 
V mesons are more than twice as heavy as the quark mass rriud = 300 - 400 McV/c^. Since the 
parameter 7 is not small, the non-relativistic reduction is not justified. Also, the very strong S-meson 
central attraction is just a replacement of the real processes of the 27r exchange, the 7rp exchange, the 
A excitations and so forth. The V mesons are supposed to behave as composite particles of the (qq) 
pairs. Furthermore, the choice of terms in Eq. (2.34) is quite ad hoc and phcnomenological. We should 
consider Eq. (2.34) an effective interaction to simulate those residual interactions between the quarks 
which are not taken into account by the FB interaction. The full expression of the Born kernel Eq. 
(2.31) is lengthy, but the direct term from the EMEP's is simple enough to show the main characteristics 
of the EMEP's introduced in the present model: 



3 



9 



An 



Ubkf 



X, 



c 

0D4 



-1 + 



2(3m„d)^ 



9 

262 



2(3m„d) 



in ■ S 



ME' 



2(3m„d) 



in 



An 



mz 



An 



k +m 



1 



(cTi ■ k){a2 ■ fe) - (1 - cs){m^ + k^)-{cri ■ (T2 



c 

0D4 



1 + 



^ (3m„,m)2^°^+ 

3 

X fc] • [0-2 X k\ 



2(3m„<i)2 
3 



9 

262 



(o-i • n)(<r2 • n) - (1 - c^.J ( — + "p" ) • ^^2) 



k' 



- 2rf 



Sniudm 



X^Ein 



-X, 



0D+ 



m 



Si-) 



(2.37) 



Here X^j-,^ represents the spin-flavor factors related to the spin-flavor operators in Eq. (2.34). Besides 
the Gaussian factor exp{— (6A;)2/3} related to the (Os)^-cluster wave functions, the finite size effect of 
the baryons also appears as the constant zero-point oscillation terms accompanied with the terms, 
which appear in the S- and V-meson contributions. For the QLS force, the same effect appears as the 
tensor force having the form [cti x k] ■ [<T2 x k] . The magnitude of this term is about one third, compared 
to the original tensor term appearing at the level of the qq interaction. The advantage of using the 
QLS force in Eq. (2.34), instead of the tensor force, is that we can avoid the n-p cancellation of the 
tensor force for the coupling term of the S and D waves. The £1 parameter of the NN interaction is 
very sensitive to this coupling strength. 

We have four QM parameters; the h.o. size parameter b for the (3g)-clusters, the up-down quark 
mass rriud, the strength of the quark-gluon coupling constant as, and the mass ratio of the strange to 



14 



up-down quarks A = {iris/mud)- A reasonable range of the values for these parameters in the present 
framework is 6 = 0.5 - 0.6 fm, m„d = 300 - 400 MeV/c^, as ~ 2, and A = 1.2 - 1.7. Note that we 
are dealing with the constituent QM with explicit mesonic degree of freedom. The size of the system 
determined from the (85) wave function with h (the rms radius of the (3g) system is equal to h) is 
related to the quark distribution, which determines the range in which the effect of the FB interaction 
plays an essential role through the quark-exchange kernel. The internal energies of the clusters should 
be calculated from the same Hamiltonian as used in the two-baryon system, and contain not only the 
quark contribution but also various EMEP contributions. The value ol as is naturally correlated with 
6, niud-i and other EMEP parameters. This implies that as is merely a parameter in our framework, 
and has very little to do with the real quark-gluon coupling constant of QCD. The quark and EMEP 
contributions to the baryon mass differences between N and A and between A and E, calculated in the 
isospin basis, is tabulated in Table III of Ref. [42]. 

For the EMEP part, we have three parameters /i, /s, and 9 for each of the S, PS, Ve (vector-electric) 
and Vm (vector-magnetic) terms. The F/ (F + D) ratio is no longer a free parameter in the present 
framework, but takes the pure SUq values «e = 1 (for S and Ve) and a.^ = 2/5 (for PS and Vm). It is 
convenient to use the coupling constants at the baryon level, in order to compare our result with the 
predictions by other OBEP models. These are related to the coupling constants at the quark level used 
in Eqs. (2.34) and (2.37) through a simple relationship 

Ji — 3fi'i , Js — 98 , h — h ^ h — 3/8 ) 

fl^'-m , fl'-fs , f^'^ = fr , fi'^^lfs'- (2.38) 

Through this replacement, the leading term for each meson in Eq. (2.37) precisely coincides with that 
of the OBEP with Gaussian form factors. In the present framework, the S-meson masses for e, S*, S 
and K, arc also considered free parameters within some appropriate ranges, expected from some real 
mesons, /o(600) or a, /o(980), ao(980) and Xq(800) or k, respectively. We further introduce three 
parameters to improve the fit of the A^A^ phase shifts to the empirical data [42] , cs the strength factor 
for the delta-function type spin-spin contact term of the one-pion exchange potential (OPEP), Cggg the 
strength factor for the spin-spin term of the QLS force, and CgT the strength factor for the tensor term 
of the FB interaction. 

We determine these parameters so as to fit the phase shift analysis SP99 [61] for the np scattering 
with the partial waves J < 2 and the incident energies Tiab < 350 MeV, under the constraint that the 
deuteron binding energy and the ^5*0 A"A" scattering length are reproduced. Some parameters are very 
sensitive to reproduce the available data for the low-energy YN total cross sections. In the parameter 
search, we minimize the value defined as 

V^=|^E(^^^^^'-^T^)'}' ' (2.39) 

where no experimental error bars are employed because they are not given in the energy-dependent 
solution of the phase-shift analysis. Here the sum over i — 1 - N covers various angular momenta and 
energies, and the mixing parameters, ei and 62, in units of degree. The value a/x^ gives a measure 
for the deviation of the calculated phase shifts from the empirical values. The result of the parameter 
search is shown in Table 2. The parameters of the previous model FSS are also shown for comparison. 
The resultant values are a/x^ ~ 3° in FSS and — 0.59° in fss2 for the np scattering, when the 
energies Tiab = 25, 50, 100, 150, 200, 300 MeV and the partial waves up to J = 2 are employed. 
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Table 2: Quark-model parameters, SU-i, parameters of the EMEP's, S-meson masses, and some reduction 
factors cs etc. for the models fss2 and FSS. The p meson in fss2 is treated in the two-pole approximation, 
for which mi and m2 (/92) are shown below the table. 





h (fm) 


mud (MeV/c2) 


as 


A = ms/mud 


fss2 


0.5562 


400 


1.9759 


1.5512 


FSS 


0.616 


360 


2.1742 


1.526 




f! 


/I 




el 1) 


fss2 


3.48002 


0.94459 


33.3295° 


55.826° 


FSS 


2.89138 


1.07509 


27.78° 


65° 




Jl 








fss2 


— 


0.26748 


— 


(no r], f]') 


FSS 


0.21426 


0.26994 


-23° 






£ VC 
Jl 


£ ve 

/s 


Jl 


Js 


fss2 


1.050 





1.000 


2.577 


(MoY/c^) 






/"<) 


/;/,, 


fss2 


800 


1250 


846^) 


936 


FSS 


800 


1250 


970 


1145 




cs 


^qss 






fss2 


0.4756^) 


0.6352 


3.139 




FSS 


0.381 









1) 0f is used only for ^N{I = 3/2). 

2) ^ 35.264° (ideal mixing) and two-pole p meson with mi (/3i) = 664.56 MeV/c^ (0.34687) and 
m2 {P2) = 912.772 MeV/c^ (0.48747) [62] are used. 

3) For the NN system, ms = 720 MeV/c^ is used. 

4) Only for tt, otherwise 1. 

5) The enhancement factor for the Fermi-Breit tensor term. 

2.3 The Lippmann-Schwinger formalism for {Sq)-{3q) RGM and the G- 
matrix equation 

As mentioned in the Introduction, the RGM equation (2.15) is solved in the momentum representation. 
The basic LS-RGM equation is given by 

r^a(p, q- E) = V^^{p, q;E) + J2 ^ / dk V^f,{p, k; £;)^^^--l-- T^„(fc, q; E) , (2.40) 

where the quasipotential V^/3(p, q; E) is calculated from the basic Born kernel (2.31) through 

V,^{k, k'; E) = ^[ M^^ik, k'- E) + {-If-Va M^^ik, -k' ■ E) ] . (2.41) 
Using the partial wave decomposition of the LS-RGM equation [44] , we introduce the definition 

T,^{k,k';E)= J2 T%^^sAk^';E)J2{^'m'S'S',\JM){emSS,\JM)Yern'(^^^^ (2.42) 

JMie' mm' 
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where 7 = [1/2(11) Ci, l/2(ll)c2] S' S'^YIIz;V' . The prime on the summation symbol indicates that i 
and £' are hmited to such values that satisfy the condition (—1)^ = {—1)^V and (—1)^ = (— 1)*^ V. We 
write the spin quantum numbers S and S' of the T-matrix explicitly for the later convenience, although 
these are already included in the definition of a and 7. The LS-RGM equation is reduced to 

f3S"i" ^ ' •^'^ 

"^prrpqr^ '^pS"e",<yseik,Q;E) , (2.43) 



X 

where E = E^"^ + (/i^/2/i^)A;| and V^g,(^,^^gf^{p,q; E) is the partial-wave decomposition of V^a{p,Q', E). 
We use the Gauss-Legendre 20-point quadrature formula to carry out the numerical integration. The 
Lippmann-Schwinger equation (2.43) involves a pole &X, k — kp in the Green function. A proper 
treatment of such a singularity for positive energies is well known. Here we use the technique developed 
by Noycs [53] and Kowalski [54] , and separate the momentum region of k (and also p and q) into two 
parts. After eliminating the singularity, we carry out the integral over Q < k < kp hy the Gauss- 
Legendre 15-point quadrature formula, and the integral over kp < k < 00 hy using the Gauss-Legendre 
20-point quadrature formula through the mapping, k = kp + tan [(7r/4)(l + x)]. 

Once we obtain the phase shift parameters for all the partial waves, it is straightforward to calculate 
the scattering amplitude, which is defined by the on-shell T-matrix through 

M,«(q^, q,- E) = - ^^J^^f ^^ T^.(g/, E) . (2.44) 

Here and q^ arc related to the total energy E by the on-shell condition E = E™^ + = E™^ + 
with = {h'^ /2fij)qf'^ and Sa = {h^ /2^ct)qi^ ■ For the calculation of the polarization obscrvables, it 
is customary to write the scattering amplitude Eq. (2.44) in the invariant form with respect to the 
space-spin operators. In the YN scattering, the three basic vectors, k — q^ — q^, q = (q^ + q^)/2 and 
n — [q^^ X qj:] — [q X k] are not mutually orthogonal. We therefore use k — q^ — q^^, n — [q x k] and 
P — [k X n] — k^q — {k ■ q)k to define the invariant amplitudes: 

M{qf, g.; E) =go + ho i{cri + cr2)-n + h_ i(<Ti - 0-2) • n /i„ (<Ti • n) • (0-2 • n) 

+ hk (CTI ■ k) ■ (0-2 ■k) + hp (<7i ■ P) ■ (0-2 ■ P) + /+ {(0-1 • fc)(0-2 ■ P) + (0-1 ■ P)(<T2 • fc)} 

+/_ { (o-i • k) {(T2 ■ P) - (o-i • P) {a2 -k)] . (2.45) 

The eight invariant amplitudes, go, - ■ ■ , /-, are complex functions of the total energy E and the scatter- 
ing angle, cos^ = (Qi-Qf)- Three of the eight invariant amplitudes, /+ and /_, do not appear for the 
A^A^ scattering due to the identity of two particles (/?._ and /__) and the time-reversal invariance (/_ and 
/+). These three terms correspond to the non-central forces characteristic in the YN scattering; i.e., h- 
corresponds to LS^~\ /+ to Si2{r,p), and /_ to LS^~^a, respectively [59]. In particular, the antisym- 
metric LS interactions, LS'^'^ and LS^~^(j, involve the spin change between and 1, together with the 
transition of the fiavor-exchange symmetry V 7^ V. In the A^A" scattering this process is not allowed, 
since the flavor-exchange symmetry is uniquely specifled by the conserved isospin: V = (— 1)^""^. On the 
contrary, these interactions arc in general all possible in the YN scattering, which gives an intriguing 
interplay of non-central forces. 

The invariant amplitudes are expressed by the yS-matrix elements through the partial-wave decom- 
position of the scattering amplitude Eq. (2.44). Expressing the invariant amplitudes as 

M^^iqf, qf, E)= R's'e,ase J2{i'm'S'S',\JM){emSS,\JM) Y,m'{qf) Y;jq,) , (2.46) 

JMW mm' 
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we obtain the partial- wave component R^s'e',ase — {^/'^'^){^^S't' ,aSi ~ ^'y,aSs',sSt',£) as follows: 

The formulas given in Appendix D of Ref. [44] are then used to reconstruct the invariant amplitudes 
by the solution of the LS-RGM equation (2.43). All the scattering observables are expressed in terms 
of these invariant amplitudes. For example, the differential cross section and the polarization of the 
scattered particle are given by 

^ = ao(^) = \go? + \ho + /i-l' + \ho -h.\^ + \K\^ + \hk\^ + \hp\^ + I/+ + /_|^ + I/+ -f.\\ 

P{e) = 2ImM/io + /i-)* + M^o-/i-)1+2Im[/ife(/+-/_)*-/ip(/+ + /_)*] . (2.48) 

The LS-RGM equation (2.40) is readily extended to the G-matrix equation by a replacement of the 
free propagator with the ratio of the angle-averaged Pauli operator and the energy denominator: 

G,^{p, q- K, u) = V,^(p, q;E) + Y,j^,jdk V,p{p, fe; E)-^^^^ Gp^{k, q; K, u) . (2.49) 

A detailed description of this formalism is given in Ref. [45] . The formula to calculate the Scheerbaum 
factor [63] for the s.p. spin-orbit potential from the G-matrix solution is given in Ref. [46]. The energy 
dependence of the quasipotential V^aiP-, q', E) in Eq. (2.49) is dealt with as follows. The total energy of 
the two interacting particles in the nuclear medium is not conserved. Since we only need the diagonal 
G-matrices to calculate s.p. potentials and the nuclear-matter properties in the lowest-order Brueckner 
theory, we simply use = E'^^ - E'^^ + {h^/2iXa)q^ both in Kya(p, Q] E) and Ky/3(p, fc; E) in Eq. (2.49). 
The meaning and the adequacy of this procedure are discussed in Ref. [55] in a simple model. 

Calculations in the isospin basis is sufficient to obtain the main characteristics of the baryon-baryon 
interactions. For more detailed investigations of the CSB etc., we also solve the LS-RGM equation 
(2.43) in the particle basis to incorporate the pion-Coulomb correction. We use the empirical baryon 
masses and evaluate spin-flavor factors for the charged pions and the neutral pion separately. The other 
spin-flavor factors for heavier mesons and the FB interaction are generated from those in the isospin 
basis through the simple isospin relations. The Coulomb force is introduced at the quark level by using 
the quark charges. The exchange Coulomb kernel has the same striicturc as the color- Coulombic term 
of the FB interaction. All the necessary expressions for the spatial integrals arc given in Ref. [42]. 

The standard technique by Vincent and Phatak [64] is employed to solve the Lippmann-Schwinger 
equation with the Coulomb force in the momentum representation. This technique requires the in- 
troduction of a cut-off radius Rc for the Coulomb interaction. In the RGM formalism, we have to 
introduce this cut-off at the quark level, in order to avoid the violation of the Pauli principle. The 
two-body Coulomb force assumed in the present calculation is therefore written as 

U?^"- = QiQje^—e{Rc - Tij) , (2.50) 

fij 

where Qi is the charge of the quark i in units of e, and G the Heaviside step function. The Coulomb 
contribution to the internal energies becomes zero for the positive charged particles like the proton and 
E+. The direct Coulomb potential is given by 

VD(i?) = ^1^26^^ jerf (/?i?) - ^ [erf (/5(i? + Rc)) + erf (/5(i? - Rc))]\ , (2.51) 
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where Z^e {i = 1, 2) is the charge of each baryon, P = \^/{2b), and erf (x) = (2/^^) e~^^dt stands 
for the error function. The exchange Coulomb kernel is also slightly modified from the pure Coulomb 
kernel. Note that, even in the np and nn systems, we have small contributions from the Coulomb 
interaction through the exchange Coulomb kernel. The value Rc should be sufficiently large to be free 
from any nuclear effect beyond Rc- Then the final S'-matrix is calculated from the condition that the 
wave function obtained by solving the Lippmann-Schwinger equation with the cut-off Coulomb force is 
smoothly connected to the asymptotic Coulomb wave function. In principle, this procedure is correct 
only when the direct Coulomb potential has a sharp cut-off. We have developed a new method [65] 
which can be applied even when the cut-off is not sharp like Eq. (2.51). In this method, we introduce 
two extra radii and Rout, which satisfy i?in < Rc < Rout- We choose i?out such that the direct 
Coulomb potential Eq. (2.51) becomes negligibly small for r > -Rout- On the other hand, R^^ is selected 
such that the Coulomb wave functions are approximately local solutions of Eq. (2.51), avoiding the 
effects of the nuclear force and the smooth transition of the cut-off effect. We solve a local-potential 
problem of Eq. (2.51) from R^ut to i?in, and prepare two independent modified plane waves, ui{k, i?in) 
and ve{k,Rin)- The standard procedure by Vincent and Phatak can now be applied to the S'-matrix 
of the LS-RGM equation by using these ^^(A;, i?in), V£{k,Rii^), in order to obtain the nuclear S'-matrix 
with the Coulomb effect. We find that this procedure gives quite accurate result within the difference 
of less than 0.1 degree, unless the incident energies are extremely small. For the practical calculations, 
we have chosen Ri^ - Rc - -Rout =6 - 8 - 12 fm. 

For the YN interaction, more consideration is required for the treatment of the threshold energies. 
Here we use the non-relativistic expressions for the threshold energies in the calculations with the isospin 
basis, while in those with the particle basis the relativistic kinematics is employed. We note that the 
mass difference of E~ and E+ is about 8 MeV and is fairly large. The Ap system has the total charge 
Z — 1 and the T,~p system Z — 0. In the YN systems, the direct Coulomb term exists only in the E+p 
channel and the E~|) channel. In calculating the EMEP contribution to the masses of the E isospin 
multiplet, we find a large cancellation between the neutral and charged pion contributions. In general, 
the pion-Coulomb correction is not sufficient to reproduce the empirical mass difference of various Bq 
isospin multiplets. The disagreement of the calculated threshold energies with the empirical values is 
a common feature of any microscopic models. Fortunately, we have a nice method to remedy this flaw 
without violating the Pauli principle. As discussed in Sect. 2.1 and in Ref. [55], we only need to add a 
small correction term AG to the exchange kernel, in order to use the empirical threshold energies. The 
same technique is also applied to the reduced mass corrections. 

2.4 Three-cluster Faddeev formalism using two-cluster RGM kernels 

To apply realistic QM baryon-baryon interactions to few-baryon systems such as the hypertriton, we 
have to find a basic three-cluster equation formulated using microscopic two-cluster quark-exchange 
kernel of the RGM. This is non-trivial not only because the quark-exchange kernel is non-local and 
energy dependent, but also because RGM equations sometimes involve redundant components due 
to the effect of the antisymmetrization, the Pauli-forbidden states. Here we propose a simple three- 
cluster equation, which is similar to the three-cluster orthogonality condition model (OCM) [66], but 
employs the two-cluster RGM kernel as the interaction potential [67]. The three-cluster Pauli-allowed 
space is constructed by the orthogonahty of the total wave functions to the pairwise Pauli-forbidden 
states. Although this definition of the three-cluster Pauli-allowed space is not exactly equivalent to the 
standard definition given by the three-cluster normalization kernel, this condition of the orthogonality 
is essential to achieve the equivalence of the proposed three-cluster equation and the Faddeev equations 
which employ a singularity- free T-matrix derived from the RGM kernel (the RGM T-matrix) [67]. 

Let us first consider a simple system composed of three identical spinless clusters, interacting via 
the two-cluster RGM kernel V^^^{e) — Vd + G + eK. We assume that there is only one Pauli-forbidden 
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state \u) for each pair of the clusters, satisfying K\u) = \u). Using the equivalence of Eq. (2.19) to Eq. 
(2.20), we can separate V^'^^{e) into two distinct parts as 

V^^^{e)^V{e)+v{e) , (2.52) 

where 

V{e) = {e- Ho) - A{e - Ho)A = e\u) {u\ + AHoA - Ho , 

v{e) = AV^^^{e)A = A{Vb + G + sK) A . (2.53) 

Note that AV{e)A = and Av{e)A = v{e); namely, v{e) represents the net effect of V^^^{e) in the 
Pauli-allowed space. Using these properties, we can express Eq. (2.19) or Eq. (2.20) as 

A[s- Ho-v(e)]Ax^O . (2.54) 

The separation of V^'^^{e) in Eq. (2.52) enables us to deal with the energy dependence of the exchange 
RGM kernel independently in the Pauli-forbidden space and in the allowed space. Let us generalize Eq. 
(2.52) to 

V{u;,e) = V{u;)+v{e) . (2.55) 

We will see that the energy dependence involved in V{lj) can be eliminated by the orthogonality 
condition to the Pauh-forbidden state. Since the direct application of the T-matrix formalism to Eq. 
(2.19) leads to a singular off-shell T-matrix [55], we first consider the subsidiary equation 

[u;-Ho-V''''''{e)]x = 0, (2.56) 

with uj ^ e, and extract a singularity-free off-shell T-matrix starting from the standard T-matrix 
formulation for the "potential" V^^^ (e) . A formal solution of the T-matrix equation 

T{uj, e) = V^^'^ie) + V^^'^{e)G^o^\uj)T{uj, e) (2.57) 

with gI^\u!) = l/{u! — Hq + iQ) is given by 

T{uj,e) = f{uj,e) + {cu-Ho)\u)^—{u\{u;-Ho), 

UJ — E 

f{u;,e) = T,{u;,e)-\l + T,{LU,e)Gi-'\u;)\\u) (m| [l + G'J+)(u;)T,(u;, e)] , (2.58) 

L {u\Gi '{lJ,6)\u) L J 

where Ty{u!,s) is defined by 

T„(u;, e) = v(s) + v(s)G^o^\u;)T^(u;, s) . (2.59) 
The basic relationship, necessary for deriving the Faddeev equations of composite particles, is 

G^o^\u;)T{u;,e) = 4+^(a;)f (u;, s) + \u)^—{u\{u; - Ho) , (2.60) 

U — 6 

with 

Gi~^\uj)T{uj,e) = G'}^\uj,e)Viuj,e) - \u){u\ . (2.61) 

Here G^^\ui,e) is given by 

gW(u;,£) = GW(a;,£) - G(^\u;,e)\u) ] (u\G(^\u;,e) , (2.62) 
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with G\;\oj, e) = 1/{uj-Hq- v{e) + iO). The full T-matrix T(w, e) in Eq. (2.58) is singular aX s — 
while T{uj,e) does not have such a singularity. For e ^ uj, T{u!,s) satisfies the relationship 

{u\G'^\u)T{u,e)\(f)) = {(f)\T{u,e)G'^\u)\u) = (2.63) 

for the plane wave solution |0) with the energy e, i.e., (e — i?o)|0) = 0. This is a direct consequence of 
more general relationship 

u)^0 , (2.64) 

which follows from Eq. (2.61). 

The three-body equation with the pairwise orthogonahty conditions is expressed as 

P[E-Ho- Kf^^(£a) - ^r^(£/3) - V^^'^'is^) ] = , (2.65) 

where Hq is the free three-body kinetic-energy operator and V^'^^{eo.) stands for the RGM kernel Eq. 
(2.52) for the a-pair, etc. The operator P projects on the Pauli-allowed space with [3] symmetry [66]. 
More precisely, we solve the eigenvalue problem 

J2K){ua\-^x) = X\-^x) (2.66) 

a 

in the model space with [3] symmetry, ^ [3], and define P as 

P = ^|*a)(*a| . (2.67) 

A=0 

The two-cluster relative energy in the three-cluster system is self-consistently determined through 

= {P^\ K + V^^'^ie^) \P^!) , (2.68) 

using the normalized three-cluster wave function with (P^|P^) = 1. Here is the free kinetic- 
energy operator for the a-pair. Expressing P^ = i'a + i'/s + i'-y with the Faddeev components ipa etc., 
we can derive the Faddeev equations [67] 

V'a = Gl^\E)fi'\E, e^){iPp + i^,) . (2.69) 

Note that Ta {E^Sa) is essentially the non-singular part of the two-cluster RGM T-matrix Eq. (2.58): 

fi'\E,e^)^f^iE-h^,e^) , (2.70) 

and that the solution of Eq. (2.69) automatically satisfies (wal'^a + '0/3 + '^7) = due to Eq. (2.64). 
Since we can also derive Eq. (2.65) from Eq. (2.69), these two equations are completely equivalent. 

This three-cluster Faddeev formalism is applied in Ref. [67] to a toy model composed of three di- 
neutrons {3d') and a more realistic model composed of three a-clusters. The di-neutron clusters are 
described with (Os)^ h.o. wave functions with z/ = 0.12 fm~^ and a Serber-type Gaussian interaction 
is assumed between two neutrons. The Pauli-forbidden state between the two di-neutrons is an A?" = 
(Os) state, where is the total h.o. quanta for the relative motion. We adjusted the strength of the 
two neutron interaction such that 2d' has a zero-energy bound state, and obtained a 3d' bound state 
of about Esd' ~ —0.5 MeV. Using the self-consistent procedure for the energy parameter e in the d'd' 
RGM kernel, we have obtained a complete agreement between this Faddeev calculation and a variational 
calculation in terms of the translationally invariant h.o. basis. 



{u\ l + G^o^\ou)T{u,e) 



= 



l + T{u;,e)Gi-^\u;) 
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In the 3a system, the structure of the 2a normahzation kernel K is not that simple. The eigenvalue 
7iv for K is given by 77V = 2^^^ — 36j^^o for the even N. In the relative 5*- wave we have two Pauli- 
forbidden states, (Os) and (Is), while in the D-wme only one (Od) h.o. state is forbidden. The relative 
motion between the two a-clusters starts from N — A h.o. quanta. The large value 74 = 1/4 in the 
Pauli-allowed space makes the self-consistent procedure through Eq. (2.68) very important. Because of 
this rather involved structure of the 2a Pauli-forbidden states, the Faddeev equations of the 3a system 
allow the existence of the redundant or trivial solutions in the Faddeev equation (2.69) [68]. These 
redundant solutions are generated from solving the overlap eigenvalue equation for the rearrangement 

{u\S\un = r\n , (2.71) 

where S is the permutation operator S — (123) + (123)^ and = is a product of two 

functions corresponding to the two momentum Jacobi-coordinate vectors, k and q, respectively. The 
eigenstate with r = — 1 gives a Faddeev component ipQ = Goluj"^), which gives a vanishing total wave 
function = (1 + S)iIj'^ = but satisfies the Faddeev equation (2.69) owing to the orthogonality 
relations of the RGM T-matrix Eq. (2.64). In the 3a system with the total angular momentum L = 0, 
we have two such solutions with t — —1. In spite of this complexity, we can eliminate these redundant 
components by slightly modifying the Faddeev equation in Eq. (2.69). For details, Ref. [68] should be 
referred to. Comparison between the Faddeev calculation and the variational calculation in the h.o. 
basis was carried out in the same condition, incorporating the cut-off Coulomb force with Rc = 10 fm. 
For the two-body effective interaction, we use the three-range Minnesota force [69] with the Majorana 
exchange mixture u — 0.94687. The h.o. width parameter, u — 0.257 fm~^, is assumed for the (Os)^ 
a-clusters. The aa phase shifts are nicely reproduced in the aa RGM calculation, using this effective 
A^A^ force. We have obtained E^a = -9.592 MeV, e = 13.481 McV, and C(o4) = 0.971, by taking the 
partial waves up to A = 8 in the Faddeev calculation [70]. Here C(o4) is the amplitude of the lowest 
shell-model component with the SU3 (04) representation. The corresponding variational calculation 
using up to = 60 h.o. quanta yields E^^ = -9.594 MeV, e = 13.480 MeV, and C(o4) = 0.971, which 
implies that the difference is only 1-2 keV. We have also compared the present Faddeev calculations 
of the 3a system with the fully microscopic 3a RGM or GCM calculations, using various A'^A'^ effective 
forces [71]. We find that the present 3a ground-state energies are only 1.5 - 1.8 MeV higher than those 
of the microscopic calculations. This implies that the genuine three-cluster exchange effect, which is 
neglected in the present three-cluster formalism, is attractive in nature, and is not repulsive claimed in 
many semi-microscopic 3a models. Oryu et al. carried out 3a Faddeev calculation using 2a RGM kernel 
[72]. They obtained very large binding energy, —E^^ = 28.2 MeV with no Coulomb force. Since the 
effect of the Coulomb force is at most 6 MeV, this 3a binding energy is extremely large. This may result 
from that they did not treat the eK term in the RGM kernel properly, and that the Pauli-forbidden 
components are not exactly eliminated. 

The present three-cluster Faddeev formalism is also applied to the three-cluster OCM in Ref. [73]. 
For this apphcation, we only need to replace V^'^^{e) in Eq. (2.52) with 

yOCM(^) ^^\u){u\+A {ho + V'''') A -ho, (2.72) 

where ho is the free two-cluster kinetic-energy operator and V^^ is an appropriate effective potential 
which reproduces the properties of the two-cluster system. If the Pauli-forbidden state \u) is a real 
bound state of V^^, satisfying (e^ — ho — V^^) \ub) = 0, Eq. (2.72) is further simplified to 

yOCM(^) ^ yeff ^ _ es)\uB){uB\ . (2.73) 

In these cases, the corresponding OCM T-matrix T{u,e) defined in Eq. (2.58) becomes £- independent. 
The Faddeev equation takes the same form as Eq. (2.69), but no self-consistency condition hke Eq. 
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(2.68) is necessary. By this procedure, we can prove that the solutions of the pairwise three-cluster OCM 
equation (2.65) (with V^'^^{e) replaced by V^^) are achieved by Kukulin's method of orthogonalizing 
pseudopotentials [74] . Here again, some complexity arises for the 3a system, when the deep aa potential 
by Buck, Priedrich and Wheatly [75] is employed with the bound-state solutions for the Pauli-forbidden 
states. In this case, one of the redundant solutions with the dominant [21] N = 4 (40) component 
does not become completely redundant, but becomes an almost redundant Faddeev component with 
the eigenvalue r ~ —0.999 in Eq. (2.71). Such a solution allows a large admixture of the [3]-symmetric 
So; component given by 

which has a large overlap with the shell- model hke compact configuration, |[3] = 8 (04)). The 
ground-state solution of the modified Faddeev equations or similar equations (see, for example, Eq. (22) 
of Ref. [76]) gives a large binding energy with this dominant configuration, and a small admixture of 
the redundant components will never be completely eliminated. On the other hand, if the Kukulin's 
method of orthogonalizing pseudopotentials is used to eliminate the redundant components completely, 
one would also repel such a component as 0r ' , resulting in a very small binding energy of the 3a ground 
state [77, 78, 79]. Refer to Ref. [76] for details. 



3 Results and discussions 
3.1 Two-nucleon system 

As already discussed, the reliability of the YN and YY interaction model is judged from whether or 
not the A^A^ data are reproduced in such a common framework that can be extended to the YN and 
YY systems. A motivation to upgrade the model FSS [40] to fss2 [42] is in a better description of the 
A^A^ data by including the important effects of the V mesons and the momentum-dependent higher- 
order terms of the S- and V-meson central forces. For example, the phase shift in FSS is more 
attractive than experiment by 10° around Tiab = 300 MeV, because the one-pion tensor force is too 
strong. The strong one-pion tensor force is partially weakened in the standard OBEP's by the p-mcson 
tensor force. We use the QLS force mainly from the p-meson exchange as a phcnomenological ingredient 
of this cancellation. Furthermore, some phase shifts of other partial waves deviate from the empirical 
ones by a couple of degrees. Another improvement is required for the central attraction. The G-matrix 
calculation using the quark-exchange kernel explicitly [45] shows that the energy-independent attraction, 
dominated by the e-meson exchange, is unrealistic, since in our previous models the s.p. potentials in 
symmetric nuclear matter show a strongly attractive behavior in the momentum region 5-20 fm~^. 
This fiaw is removed in Ref. [44] by introducing the momentum-dependent higher-order term of the 
fiavor-singlet S-meson exchange potential, which was first carefully examined by Bryan and Scott [80]. 
In the higher energy region, the LS term of the S mesons also makes an appreciable contribution. 

Figures 1(a) - 1(h) compare the np phase shifts and the mixing angles ej predicted by fss2 with those 
of SP99 by Arndt et al. [61]. For comparison, the FSS results are also shown with the dotted curves. 
The phase-shift parameters are shown for the partial waves up to J = 4 with the energies less than 350 
MeV. More information is obtained from the QMPACK homepage [81]. For the non-relativistic energies 
less than Tiab = 300 MeV, the A^A^ phase-shift parameters are almost perfectly reproduced, except for 
the state. The values in Eq. (2.39) are calculated using the phase-shift analysis SP99 [61] 

and PWA93 [82] of the Nijmegen group. The fss2 gives = 0.59° and 0.60° for SP99 and PWA93, 
respectively. These are compared to the values of CD-Bonn [47] and ESC04 [26]; a/x^ = 0.52° 
(0.16°) and 1.16° (1.01°) for CD-Bonn and ESC04, respectively, when SP99 (PWA93) is used. 
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Figure 1: Calculated np phase shifts by fss2 in the isospin basis, compared with the phase-shift analysis 
SP99 by Arndt et al. [61] The dotted curves indicate the results given by FSS. 
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Figure 2: Decomposition of the ^Pj phase shifts 
for the np scattering to the central (Sc), LS (Sls) 
and tensor {6t) components up to the energy 
Tiab = 800 MeV. The results given by fss2 (sohd 
curves) and FSS (dashed curves) are compared 
with the decomposition of the empirical phase 
shifts SP99 [61]. 
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For energies higher than 300 MeV, the inelastic process of the pion production takes place, and the 
effect of the AN channel becomes important. The inelasticity parameters determined from the phase 
shift analysis, however, show that they are still rather small even up to the energies Tiab = 800 MeV. We 
have therefore examined the results of the single-channel NN calculation up to 800 MeV in Ref . [42] , in 
order to investigate a possible effect of the AN channel coupling. Up to this energy, the phase shift 
is well reproduced by the QLS component. Even in the other partial waves, the improvement of the 
phase-shift parameters is usually achieved. This includes 1) ^Pq, ^-Pi, and ^G^ phase shifts, 2) ^5"!, ^-S'o, 
^Pi, ^Fs, and ^H^ phase shifts at higher energies Tiab = 400 - 800 MeV, and 3) some improvement in 
phase shift and £2 mixing parameter. On the other hand, '^P2 and ^D^ phase shifts turn out worse 
and ^F^ phase shift is not much improved. The problem of fss2 lies in the ^P2 and ^D^ phase shifts 
at the intermediate and higher energies Ti^b = 300 - 800 MeV. The empirical ^P2 phase shift gradually 
decreases if we ignore the weak dispersion-like behavior. Our result, however, decreases too rapidly. 
Our "^Do, phase shift is too attractive by 4° - 6°. The disagreement of the '^Ds phase shift and the 
deviation of the ^Di phase shift at higher energies suggests that our description of the central, tensor 
and LS forces in the "^E states requires further improvement. The insufficiency in the ^'O partial waves 
is probably related to the imbalance of the central force and the LS force in the short-range region. 
The decomposition of the ^Pj phase shifts to the central, LS and tensor components, shown in Fig. 2, 
implies that the central force is too repulsive at higher energies Tiab > 400 - 500 MeV. We find that 
whenever the discrepancy of the phase shifts between the calculation and the experiment is large, the 
inelasticity parameters are also very large. In particular, the inelasticity parameters of the ^P2, ^D2, 
and states rise very rapidly as the energy increases, and reach more than 20° at Tiab = 800 MeV. 
The elastic phase shift for each of these states shows a dispersion-hke resonance behavior in the energy 
range from 500 MeV to 800 MeV. These are the well-known di-baryon resonances directly related to 
the AA^ threshold in the isospin 1 = 1 channel. The present single-channel calculation is not capable 
of describing these resonances, unlike the ^ So{N N)-^ Do{N A) CCRGM calculation in Ref. [83]. 

Figure 3 compares with experiment [61] the fss2 (solid curves) and FSS (dashed curves) predictions 
of the differential cross sections [da/dVt) and the polarizations {P{9)) for the elastic np scattering. 
The same observables for the elastic pp scattering are plotted in Fig. 4, and some other observables 
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np differential cross sections and the polarization compared with the experiment 
is performed using fss2 (solid curves) and FSS (dashed curves) with the full pion- 
included. 
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ition {D{9)), the vector analyzing powers (^(6'), A'{9), R{0)), etc., in Figs. 5 and 6. 
are obtained in the particle basis with the full Coulomb force incorporated. We find 
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Figure 5: The same as Fig. 3 but for the other spin observables. 



some improvement in the differential cross sections. First, the previous FSS overestimation of the np 
differential cross sections at the forward angle at Tlab > 100 MeV is corrected. Secondly, the bump 
structure of the np differential cross sections around ^cm = 130° at energies Tiab = 300 - 800 MeV has 
disappeared. The overestimation of the pp differential cross sections at 9cm = 10° - 30° at energies 
Tiab = 140 - 400 MeV is improved. However, the essential difficulties of FSS, namely the oscillatory 
behavior of the np polarization around O^m = 110° and that of the pp polarization around the symmetric 
angle 6cm ~ 90° for higher energies Tiab > 400 MeV are not resolved. Furthermore, the pp differential 
cross sections show a deep dip at angles 9cm < 30° and > 150° for Tiab > 500 MeV. The low-energy pp 
cross sections at = 90° for Tiab ^ 100 MeV are still overestimated. 

In order to find a possible reason for the unfavorable oscillations of our polarizations, we compared 
in Fig. 3 of Ref. [42] the five independent pp invariant amplitudes at the highest energy Tiab = 800 MeV, 
predicted by fss2 and the phase-shift analysis SP99. They are composed of the real and imaginary parts 
of Qq (spin-independent central), /iq {LS), hn ( («Ti • n)(«T2 • n)-type tensor), ( («Ti • k){cr2 ■ k)- 
typc tensor), and hp ( (cti ■ P)(cr2 ■ P)-type tensor) invariant amplitudes in Eq. (2.45). The result 
by SP99 was calculated using only the real parts of the empirical phase-shift parameters. Recalling 
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Figure 7: The momentum dependence of the s.p. potentials Ub{(1) predicted by the G-matrix calculation 
of fss2 in the continuous prescription for intermediate spectra: (a) The real part, Re Ub{(1)- (b) The 
imaginary part, Im Ub{(i)- The nucleon s.p. potentials obtained by the f^^ p prescription are also shown 
for fss2, the Paris potential [84] and the empirical phase shifts SP99 [61]. The momentum points selected 
correspond to Tlab = 100, 200, 400, 800, and 1600 MeV for the NN scattering. The partial waves up 
to J = 8 are included in fss2 and the Paris potential, and J = 7 in SP99. 



that the polarization is given by the cross term contribution of the central, LS, and tensor invariant 
amplitudes (i.e., P{e) = 2 lm[{go + hr,){ho)*] from Eq. (2.48) since /i_ = /± = for A^A^), we find 
that the disagreement in Im and Re Hq is most serious. Since the oscillatory behavior of Im hn in 
SP99 also appears in Im hk and Im hp, it is possible that this is an oscillation caused by the NN-AN 
channel coupling through the one-pion spin- spin and tensor forces. Prom the figure, we can also expect 
the underestimation of the differential cross sections at ^cm < 30° in this energy region. Namely, the 
imaginary part of 5^0 is too small and the real part of qq is strongly reduced in fss2. 

Another application of the invariant amplitudes is the f^^p prescription for calculating the s.p. poten- 
tials of the nucleons and hyperons in nuclear matter. It is discussed in Ref. [44] that the s.p. potentials 
predicted by the model FSS in the G-matrix calculation show fairly strong attractive behavior in the mo- 
mentum interval g = 5-20 fm~^ for all the baryons. In particular, Uf^{q) in the continuous prescription 
becomes almost —80 MeV at g = 10 fm~^. This momentum interval corresponds to the incident- energy 
range Tiab = 500 MeV -8 GeV in the A^A^ scattering. The f^^p prescription is a convenient way to eval- 
uate the s.p. potentials in the asymptotic momentum region in terms of the spin-independent invariant 
amplitude at the forward angle go{9 = 0). Since the present model fss2 incorporates the momentum- 
dependent Bryan-Scott term, the asymptotic behavior of the s.p. potentials in the large momentum 
region is improved. We can see this in Fig. 7, where the s.p. potentials of A^, A, and S calculated in 
the G-matrix approach are shown in the momentum range q = - 10 im-\ Fi gures 7(a) and 7(b) show 
the real and imaginary parts of UB{q), obtained in the continuous choice for intermediate spectra. The 
solid curves for the nucleon s.p. potential are compared with the results by the f^^p prescription with 
respect to the T-matrices of fss2, the Paris potential [84], and the empirical phase shifts SP99 [61]. 
The partial waves up to J < 8 are included in fss2 and the Paris potential, and J < 7 in SP99. The 
momentum points calculated correspond to the energies Tiab = 100, 200, 400, 800, and 1,600 MeV. We 
find that the real part of UN{q) nicely reproduces the result of the G-matrix calculation even at such a 
low energy as Tiab = 100 MeV. On the other hand, the imaginary part by the t^^p prescription usually 
overestimates the exact result especially at lower energies. 
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The deuteron properties are calculated by solving the LS-RGM equation with respect to the relative 
wave functions /o(/c) and f2{k) in the momentum representation (see Appendix B of Ref. [42]). The 
properly normalized wave functions in the Schodinger picture are not fi{k) but — \fN where 
N represents the normalization kernel [40]. The (S-wave and D-wave wave functions in the coordinate 
representation, u{R) and w{R), arc then obtained from the inverse Fourier transform of Fi{k). This 
process is most easily carried out by expanding F£{k) in a series of Yukawa functions ^/^j-Kkj [k^ + 7j^) 
in the momentum representation (see Appendix D of Ref. [47]). We choose 7j = 7 + (j — l)7o (j = 1 
- 11) with 7o = 0.9 fm~^. The 7 is the yS-matrix pole q — — i7, from which the deuteron energy is 
most accurately calculated by using the relativistic relation 



72 + 



(3.1) 



Figure 8 shows the deuteron wave function of fss2 in the coordinate and momentum representations, in 
comparison with that of the Bonn model-C potential [85] (dotted curves)^. We find that the difference 
between the two models is very small. Table 3 compares various deuteron properties calculated in three 
different schemes: the isospin basis and the particle basis without or with the Coulomb force. They are 
also compared with the empirical values and the predictions by the Bonn model-C potential. The final 
value of the deuteron binding energy for fss2 is = 2.2309 MeV. If we use the non-relativistic energy 
expression,^ td = (j^/Mn) for 72 = 0.053762 fm~^ in the fuU calculation, we obtain = 2.2295 MeV 
and the difference is 1.4 keV. The differences within the deuteron parameters calculated in the three 
different schemes are very small, except for the binding energy 64. In particular, the exchange Coulomb 
kernel due to the exact antisymmetrization at the quark level gives an attractive effect to the binding 
energy, and increases by 4.8 keV. This is even larger than the relativistic correction included in Eq. 
(3.1). The deuteron D-state probability is Pd = 5.49 % in fss2, which is slightly smaller than 5.88 % 
in FSS [40]. These values are close to the Bonn model-C value Po = 5.60 % [85]. The asymptotic D/S 
state ratio r] and the rms radius arc very well reproduced. On the other hand, the quadrupole moment 
is too small by about 5 - 6%. There are some calculations [87, 88] which claim that the effect of the 
meson-exchange currents on the deuteron quadrupole moment is as large as AQ^ — 0.01 fm^. It is 
noteworthy that the Bonn model-C almost reproduces the correct quadrupole moment, in spite of the 
fact that the D-state probability is very close to ours. (On the other hand, the quadrupole moment of 



^The results of the Bonn model-C potential in Fig. 8 and in Table 3 are based on the parameterized deuteron wave 
functions given in Table C.4 of Ref. [85]. 

^In Table 3, the value of in the isospin basis is calculated using this non-relativistic formula. 



c 
.0 

"-4— • 

o 
c 

3 



> 

TO 



C 

o 



3 

Q 




c 
.0 

c 
3 



0) 

> 
CO 



c 
o 



3 

o 




Figure 8: The deuteron wave functions predicted by fss2 (solid curves) and by Bonn model-C [85] in 
the coordinate (a) and momentum (b) representations. 
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Table 3: Deuteron properties by fss2 in three different calculational schemes, compared with the pre- 
dictions of the Bonn model-C potential [85] and the experiment. 





isospin basis 


particle basis 


Bonn C 


Expt. 


Ref. 






Coulomb off 


Coulomb on 








ea (MeV) 


2.2250 


2.2261 


2.2309 


fitted 


2.224644 ± 0.000046 


[86] 


Pd (%) 


5.490 


5.490 


5.494 


5.60 






r] = Ad/ As 


0.02527 


0.02527 


0.02531 


0.0266 


0.0256 ± 0.0004 


[89] 


rms (fm) 


1.9598 


1.9599 


1.9582 


1.968 


1.9635 ± 0.0046 


[86] 


Qd (fm2) 


0.2696 


0.2696 


0.2694 


0.2814 


0.2860 ± 0.0015 


[90] 




0.8485 


0.8485 


0.8485 


0.8479 


0.85742 





CD-Bonn [47] is Qd — 0.270 fm^ with a smaller value Pb = 4.85%.) For the magnetic moment, precise 
comparison with the experimental value requires a careful estimation of various corrections arising from 
the meson-exchange currents and the relativistic effect, etc. 

Table 4 hsts the S- and P-wave effective range parameters for the A^A^ system, calculated in the three 
schemes. Since the pion-Coulomb correction is not sufficient to explain the full CIB effect existing in 
the np and pp ^Sq states, a simple prescription to multiply the flavor-singlet S'-meson coupling constant 
fi by a factor 0.9949 is adopted to reduce too large attraction of the pp central force. (This prescription 
is applied only to the calculation in the particle basis.) The underlined values of the scattering length 
a in Table 4 indicate that they are fitted to the experimental values. We find that the pion-Coulomb 
correction in the np ^Sq state has a rather large effect on a. The value a — —23.76 fm in the isospin 
basis changes to a = —27.39 fm due to the effect of the pion mass correction and the explicit use of 
the neutron and proton masses. It further changes to a = —27.87 fm due to the small effect of the 
exchange Coulomb kernel. These changes should, however, be carefully reexamined by readjusting the 
binding energy of the deuteron in Table 3. We did not carry this out, since the reduction of ff to 
fit these values to the empirical value a — —23.748 ± 0.010 fm does not help much to reproduce the 
CIB of the pp channel. This shortcoming might be related to the insufficient description of the low- 
energy pp differential cross sections around ^^cm = 90°, as seen in Fig. 4. It was also pointed out by 
the Nijmegen group [94] that the Coulomb phase shift should be improved by the effects of two-photon 
exchange, vacuum polarization, and magnetic moment interactions, in order to describe the ^5*0 phase 
shift precisely at energies less than 30 MeV. Neither of these effects, nor more sophisticated charge- 
symmetry breaking operators due to the interference of QED and QCD [95], are taken into account in 
the present calculation. The P-wave effective range parameters are also given in Table 4, to compare 
with a number of empirical predictions. The parameters of ^P2 state are not given, since the effective 
range expansion of this partial wave requires a correction term related to the accidental p^ low-energy 
behavior of OPEP [96]. 

3.2 Hyperon-nucleon interactions 

3.2.1 The spin-flavor SUq symmetry in the B^Bs interactions 

Our (3g)-cluster wave functions in the (3g)-(3g) RGM formalism are the strict SUq wave functions in the 
spin-flavor 5*^76 QM. This fact and the complete antisymmetrization of quarks in the RGM framework 
make it possible to investigate interesting relationship of the various BsBs interactions under some 
plausible assumptions and special roles of pions. First it should be noted that our phenomenological 
confinement potential does not have any fiavor dependence. The contribution of the r^-type confine- 
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Table 4: Effective range parameters of fss2 for the NN interactions. For the pp and nn systems, the 
calculation in the particle basis uses x 0.9949, in order to incorporate the effect of the charge- 
independence breaking (CIB). Unit of length is in fm2^+^ (a), fm-2^+^ (r) and fm"^^ (P) for the partial 
wave i. The experimental values are taken from [86], [91], [92], [93], [94], and [47]. 









isospin basis 


particle basis 


Expt. 










Coulomb off 


Coulomb on 








a 




-17.80 


-7.810 


-7.819 ±0.0026 


pp 


'So 


r 
P 




2.680 
0.061 


2.581 
-0.064 


2.794 ± 0.0014 






a 




-2.876 


-3.117 


-4.82 ±1.11, -2.71 ±0.34 


pp 




r 
P 




3.834 
-0.011 


3.735 
-0.034 


7.14 ±0.93, 3.8 ±1.1 






a 




1.821 


1.992 


1.78 ±0.10, 1.97 ±0.09 


pp 


'Pi 


r 
P 




-8.162 
0.001 


-7.5666 
0.001 


-7.85 ±0.52, -8.27 ±0.37 






a 




-18.04 


-18.04 


-18.05 ±0.3, -18.9 ±0.4 


nn 


'So 


r 
P 




2.677 
0.061 


2.677 
0.061 


2.75 ±0.11 






a 




-2.881 


-2.881 




nn 


'Po 


r 
P 




3.825 
-0.011 


3.825 
-0.011 








a 




1.823 


1.823 




nn 


'Pi 


r 
P 




-8.154 
0.001 


-8.155 
0.001 








a 


-23.76 


-27.39 


-27.87 


-23.748 ± 0.010 


np 


'So 


r 


2.588 


2.532 


2.529 


2.75 ± 0.05 






P 


0.059 


0.052 


0.052 








a 


-2.740 


-2.466 


-2.466 




np 


'Po 


r 


3.869 


3.930 


3.930 








P 


-0.013 


-0.018 


-0.018 








a 


5.399 


5.400 


5.395 


5.424 ± 0.004 


np 


'Si 


r 


1.735 


1.735 


1.734 


1.759 ±0.005 






P 


0.063 


0.063 


0.063 








a 


2.824 


2.826 


2.826 




np 


'Pi 


r 


-6.294 


-6.300 


-6.300 








P 


-0.006 


-0.006 


-0.006 








a 


1.740 


1.582 


1.582 




np 


'Pi 


r 


-8.198 


-8.185 


-8.185 








P 


0.001 


0.000 


0.000 





ment potential to the baryon-baryon interaction is zero due to the subtraction of the internal energy 
contribution. The FB interaction involves the FSB only through the quark-mass dependence, but the 
effect is rather weak. Thus we can assume that our QM Hamiltonian is approximately SU3 scalar. 
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On the other hand, the EMEP contribution is also approximately SUs scalar, except for the pion con- 
tribution. This is because the SU3 relations of the coupling constants are automatically incorporated 
in the QM description. The usage of the empirical masses for the PS and V mesons breaks the SU3 
symmetry, together with the baryon masses. The S-meson masses are varied to fit the NN and YN 
data, but these are also assumed to be SU^ symmetric in the first approximation. After all, our QM 
Hamiltonian with the EMEP's is approximately S'f/3 scalar as a whole, except for the important roles 
of pions with the very small mass. If the Hamiltonian is approximately SU3 scalar, we can expect 
that the baryon-baryon interactions with the same SU^ representation should be very similar, since 
the dependence on the internal quantum numbers becomes negligible. Prom this discussion, we can 
correlate different B^Bs interactions according to the SU3 representation of the BgBg systems. 

Table 5 illustrates how the two-baryon systems in the isospin basis, defined by Eq. (2.8), are classified 
as a combination of the SU^ bases. Entirely different structure appears between the flavor-symmetric 
and antisymmetric cases. In the ^Sq state, for example, we have several systems containing the (22) 
component dominantly. The 5'-state B^Bg interactions for these systems are expected to be similar 
to that for the A^A^ ^5*0 state. The (ll)s component is completely Pauli forbidden for the most com- 
pact (Os)^ configuration, and is characterized by the strong repulsion originating from the quark Pauli 
principle. The (00) component in the i7-particle channel is attractive because of the color-magnetic 
interaction. On the other hand, in the A"A^ sector the ^Si state with / = is composed of the pure 



Table 5: The relationship between the isospin basis and the flavor-S'C/a basis for the -BgSg systems. The 
Hmor-SU^ symmetry is given by the Elliott notation (A//). The heading V denotes the flavor-exchange 
symmetry, S the strangeness, and / the isospin. 



s 


BsBs (/) 


V — +1 (symmetric) 


V — —1 (antisymmetric) 






or ^0 


^E or ^0 





NN {I = 0) 
NN {I = 1) 


(22) 


(03) 


-1 


AA^ 
SAT (/ = 1/2) 
EAT (/ = 3/2) 


^[(11), + 3(22)] 
;^[3(ll).-(22)] 

(22) 


^[-(11), + (03)] 

^[(ll)a+(03)] 

(30) 




AA 

SA^ (I = 0) 


7l(ll)^ + iA3(22) + ^(00) 
^(11). - ^(22) + ^(00) 


(ll)a 




EN (I = 1) 


^(11)3+^1(22) 


^[-(ll)„ + (30) + (03)] 


-2 


EA 




73PO) - (03)] 




EE (J = 0) 
EE (/ = 1) 
EE (/ = 2) 


yi(ll)^-i7T3(22)-y|(00) 
(22) 


^[2(ll)„ + (30) + (03)] 


-3 


SA 

SE (J = 1/2) 
SE (/ = 3/2) 


^[(11), + 3(22)] 
;fe[3(ll).-(22)] 
(22) 


^[-(11), + (30)] 

^[(ll)a + (30)] 

(03) 


-4 


EE (/ = 0) 
EE (/ = 1) 


(22) 


(30) 
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(03) state, and the deuteron is bound in this channel owing to the strong one-pion tensor force. This 
SU3 state in the flavor-antisymmetric case is converted to the (30) state in the side of more negative 
strangeness. Since the (30) state is almost Pauh forbidden with the eigenvalue = 2/9 of the normal- 
ization kernel, the interaction is strongly repulsive. Therefore, the SS interaction is not so attractive 
as NN, since they are combinations of (22) and (30). The other SU^ state (ll)a turns out to have 
very weak interaction. After all, the strangeness S = —2 sector is most difficult, since it is a turning 
point of the strangeness. It is also interesting to see that the SE channel with 7 = 3/2 should be fairly 
attractive, since the channel contains the same (22) and (03) SU^ states as the NN system. 

Figure 9(a) shows fss2 predictions for the ^Sq phase shifts of various BgBg interactions having the 
pure (22) configuration. The corresponding phase-shift equivalent local potentials at piab = 200 McV/c, 
derived by the WKB-RGM technique [97, 98], are depicted in Fig. 10(b). Although the SS interaction 
with isospin / = 2 is very similar to the NN interaction, the other interactions generally get weaker as 
more strangeness is involved. This is a combined effect of both the FSB in the quark sector and the 
EMEP contributions. In particular, the SS interaction shows the lowest rise of the phase shift, which 
is less than 30°. Accordingly, the SS total cross sections become much smaller than the other systems, 
as seen in Fig. 9(b). This figure also shows the total cross sections of the np and S~E~ scatterings, 
in which the flavor-symmetric states have the SU^ (22) representation and the antisymmetric states 
(03) representation. The larger cross sections of the np scattering are due to the very strong effect of 
pions. There is no pion effect in the exchange Feynman diagram for the S~E~ scattering. More detailed 
analysis shows that the NN interaction is the strongest and has the largest cross sections among any 
combinations of the B^s. 

This kind of analysis is also applicable to the flavor-antisymmetric conflgurations. Figure 10(a) 
shows the phase-shift behavior of the ^Si states for the NN{I = 0) system with (03), SE(/ = 3/2) with 
(03), EN{I = 0) with (11)„, EA^(7 = 3/2) with (30), and SS(7 = 0) with (30). First we find that the 
(03) states arc attractive similarly to the (22) state. The existence of the deuteron in the NN{I = 0) 
system is due to the very strong tensor force of the one-pion exchange. Next, the (ll)a state represented 
by the S7V(/ = 0) system is weakly attractive, and the phase-shift rise is only 5°. Finally, the (30) 
states represented by the T,N(I = 3/2) and EE{I = 0) systems are both repulsive, but the repulsion of 
the SS(7 = 0) state is shghtly weaker than that of the T,N{I — 3/2) state due to the FSB. 

It should be kept in mind that the systematics based on the relationship in Tabic 5 serves as a guide 
and the deviation from the pure SU3 symmetric rules is very much model dependent. For example, 
the model FSS gives ^5*0 bound states for the EE(/ = 2), SE(/ = 3/2), and SS(/ = 1) systems with 




Figure 9: (a) ^So phase shifts for the i?8 Bg interactions with the pure (22) state, predicted by fss2. (b) 
Total elastic cross sections for the pure (22) state {pp, E^E~, S~5~) and for the (22)-|-(03) states {np, 
S~E~), predicted by fss2. The Coulomb force is neglected. 
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Figure 10: (a) ^5*1 phase shifts for some Sg Bg interactions, predicted by fss2. The SU^ configurations 
are (03) for NN{I = 0) and = 3/2), (11)„ for EN{I = 0), (30) for SiV(/ = 3/2) and EE{I = 0). 

The NN{I = 0) phase shift predicted only by the central interaction is also shown by crosses, (b) 
phase-shift equivalent local potentials of NN{I = 1), EN {I = 3/2) with Os = 33.3295° (see Table 
2), EE(7 = 2), SE(7 = 3/2), and SS(7 = 1) systems predicted by fss2 at piab = 200 MeV/c. 

the pure (22) representation. This is due to the imbalanced reduction in the short-range repulsion of 
the color-magnetic interaction and the intermediate-range attraction of the flavor-singlet scalar-meson 
exchange EMEP with increasing strangeness. Since no such bound states are observed experimentally, 
the above discussions based on the model fss2 seem to be more reliable and give an important guide 
line to construct realistic models for the BgBs interactions. 

3.2.2 S+p system 

Sine the 5*- wave phase shifts of the EN {I = 3/2) system are already given in Figs. 9(a) and 10(a), 
we only show the P-wave phase shifts in Fig. 11. The phase shifts given by FSS are also shown for 
comparison. The and states of the EN(I — 3/2) system belong to the (22) representation in the 
flavor SUs symmetry, in exactly the same as the isospin 7 = 1 A^A^ system. The phase-shift behavior of 
these partial waves therefore resembles that of the A^A^ system, as long as the effect of the FSB is not 
significant. Figure 9(a) shows that the attraction in the ^5*0 state is much weaker than that of the A^A^ 
system, and the phase-shift peak is about 26° around ps = 200 MeV/c. The ^Pj phase shifts show the 
characteristic energy dependence observed in the NN phase shifts, which is caused by different roles of 
the central, tensor and LS forces. The appreciable difference in the two models fss2 and FSS appears 
only in the state. It is discussed in Ref. [99] that the attractive behavior of the phase shift is 
closely related to the magnitude of the S+p polarization P{6) at the intermediate energies. The more 
attractive the phase-shift behavior, the larger P{9 — 90°). In Fig. 12(a), the polarization observable for 
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Figure 11: The P-wave phase shifts of the 
UN {I = 3/2) system, predicted by fss2 (sohd 
curves) and FSS (dashed curves). 
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Figure 12: S+p (a) and Ap (b) asymmetries at 
ps = 900 MeV/c, predicted by fss2 and FSS. The 
experimental data are taken from [100, 101]. 



the E+p elastic scattering at ps = 900 MeV/c is shown for the models fss2 and FSS. The experimental 
data [100] imply the asymmetry parameter a"^^"^ = 0.44 ± 0.2 at ps = 800 ± 200 MeV/c, which is not 
inconsistent with our QM predictions provided that the scattering angle is not specified. More recent 
data [101] at ps = 700 - 1, 100 MeV/c are plotted in Fig. 12(a). 

On the other hand, the and states of the T,N{I = 3/2) system have the (30) symmetry. We 
have no information on the properties of these states from the NN interaction. In our QM framework, 
however, there is very little ambiguity in the phase shifts of these states. Since the (Os)^ configuration 
in the ^'Si state is almost Pauli forbidden, the interaction in the state is strongly repulsive. On the 
other hand, the ^Pi phase shift is weakly attractive, which is caused by the exchange kinetic-energy 
kernel owing to the Pauli principle. 

The effective range parameters of the YN scattering in the single-channel analysis are given in Table 
6 with some empirical values. For the E+p system, the empirical values given in Ref. [48] should be 
compared with the results in the particle-basis calculation including the Coulomb force. We find a 
reasonable agreement both in the ^5*0 and ^5*1 states. 

3.2.3 AN system 

The total cross section for the Ap elastic scattering predicted by fss2 in the isospin basis is displayed in 
Fig. 19(e) later, together with the previous result given by FSS. The both models reproduce experimental 
data equally well in the low momentum region pa < 300 MeV/c. The total cross section has a cusp 
structure at the EA^ threshold due to the strong AN-HN{I = 1/2) ^5*1 - ^Di coupling caused by the 
tensor force of the OPEP. To see the dominant role of the tensor force of the AA^-EiV coupling in more 
detail, we show the ^5*1 and ^5*0 phase shifts in Fig. 13 calculated in the isospin basis. A cusp structure 
at the EA" threshold is apparent in the ^Si channel, while very small in the ^Sq channel. In both models. 
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Table 6: ^Sq (og, r^) and ^Si {at, rt) effective range parameters of fss2 for -Bg-Bg interactions. Unit of 
length is in fm. The experimental values are taken from [48] for E+p, and from [108] and [109] for Ap. 







isospin basis 


particle basis 


Expt. 








Coulomb off 


Coulomb on 






tts 


-2.51 


-2.48 


-2.28 


-2.42 ±0.30 


E+p 




4.92 


5.03 


4.68 


3.41 ±0.30 




at 


0.729 


0.727 


0.826 


0.709 ± 0.001 




n 


-1.22 


-1.31 


-1.52 


-0.783 ± 0.003 




as 


-2.59 


-2.59 


-2.59 


-1.8, -2.0 


Ap 


Ts 


2.83 


2.83 


2.83 


2.8, 5.0 




at 


-1.60 


-1.60 


-1.60 


-1.6, -2.2 
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3.01 


3.00 


3.00 


3.3, 3.5 


E+E+ 
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Ts 


2.34 
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2.26 
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-0.821 
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rs 
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-1.08 






Ts 
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0.262 
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2.15 


2.14 
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-4.63 


-4.70 


-3.46 
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2.39 


2.37 


2.04 






at 


-3.48 
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-2.81 
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2.52 


2.52 
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-1.43 


-1.43 






Ts 


3.20 


3.17 


3.17 






at 


3.20 


3.20 


3.20 






n 


0.218 


0.218 


0.217 






as 




-1.43 


-1.35 






Ts 




3.21 


2.70 





the attraction in the ^5*0 state is stronger than that in the '^^i state. However, the relative strengths 
of attraction in the ^Sq and ^Si states are different. The model FSS gives 5Q-Sq) — 5{^Si) ^ 30° at 

Pe ^ 200 MeV/c, while fss2 5CSo) - 5(^Si) ^ 10°. The difference of 5CSo) - 5{^Si) < 10° is required 
from many few-body calculations [102, 103, 104, 105, 106] of the s-shell A-hypernuclei. In Sect. 3.4.2, 
we will show Faddeev calculations for the hypertriton [107], using fss2, which predict that the desirable 
difference is only 0-2°. 

As seen in Fig. 19(e) later, FSS predicts especially large enhancement of the total cross sections 
around the cusp at the EA^" threshold. This is due to a rapid increase of the AA^" ^Pi - AN ^Pi transition 
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around the threshold. Figure 14 displays the S'-matrix % = r/i^e^^^'i for the AiV-SA^(/ = 1 /2) ^Pi - 
^Pi channel coupling for fss2 (upper) and FSS (lower). The FSS shows that the transmission coefficient 
r)2i, which corresponds to the AA^ ^Pi — > AA" ^Pi transition, increases very rapidly around the EA" 
threshold as the energy increases. This increase of the 7^21 (and the resultant decrease of the reflection 
coefficient riu) is a common feature of our models. The strength of the transition, however, has some 
model dependence. The transition is stronger in FSS than in fss2. The behavior of the diagonal phase 
shifts is largely affected by the strength of the ^Pi - ^Pi channel coupling (which is directly reflected 
in 7721). The model fss2 yields a broad resonance in the EA" ^Pi channel, while the FSS with stronger 
channel coupling yields no resonance in this channel. Instead, a step-like resonance appears in the AA^ 
^Pi channel. The location of the resonance is determined by the strength of the LS**^"-* force and the 
strength of the attractive central force in the SA^(/ =1/2) channel. In the QM, the central attraction 
by the S mesons is enhanced by the exchange kinetic-energy kernel, which is attractive in the SA^ ^Pi 
channel due to the effect of the Pauli principle. The LS and LS'^'^ forces also contribute to increase 
this attraction. In FSS, a single channel calculation for the UN (I — 1/2) system yields a resonance in 
the ^Pi state. When the channel coupling to the AA^ system is introduced, this resonance moves to the 
AA^ channel due to the strong LS^^^ force generated from the FB interaction. On the other hand, the 
role of the FB LS^~^ force becomes less significant in fss2, since the S- meson exchange includes only 
the LS contribution and no LS^~^ contribution. This is the reason the resonance remains in the EA^ 
^Pi channel in fss2. 

In spite of these quantitative differences in the coupling strength, the essential mechanism of the 
AN-TiN{I = 1/2) ^Pi - ^Pi channel coupling is the same for all of our models. It is induced by the 
strong FB LS^~^ force, which directly connects the two SUs configurations with the (ll)s and (ll)a 
representations in the P-wave 1—1/2 channel. In order to determine the detailed phase-shift behavior 
of each channel, including the position of the P-wave resonance, one has to know the strength of the 
LS^~'' force and the strength of the attractive central force in the SA^(/ = 1/2) channel. This is possible 
only by the careful analysis of rich experimental data for the Ap and E^p scattering observables in the 
EA^ threshold region. The asymmetry parameter data for the Ap elastic scattering at pa = 700 - 
1, 100 MeV/c are plotted in Fig. 12(b) [101]. 

Let us discuss another important feature of the AA" phase shifts induced by the P-wave coupling 
due to the LS"*^"-* force. The AA^ ^Pi and ^Pi phase shifts in Fig. 14 show weakly attractive behavior 
below the EA^ threshold. This is due to the dispersion-like (or step-like for the AA^ ^Pi state in FSS) 
resonance behavior with a large width. This attraction in the P state can be observed by the forward to 




Figure 13: The AN ^Si and ^Sq phase shifts calculated by fss2 (a) and by FSS (b). 
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Figure 14: The ^'-matrix of the AiV-EiV(/ = 1/2) ^Pi - ^Pi CC system, predicted by fss2 (upper) and 
FSS (lower) in the isospin-basis. The diagonal phase shifts Sa (left-hand side), and the reflection and 
transmission coefficients rjn (right-hand side) defined by Sij — rjije'^'^^^^ are displayed. The channels are 
specified by 1: AA^ ^Pi, 2: AA^ ^Pi, 3: EA^ ^Pi and 4: EA^ ^P^. 

backward ratio (F/B) of the Ap differential cross sections. Both fss2 and FSS give F/B > 1 below the 
EA^ threshold [41], which implies that the P-state AA^ interaction is weakly attractive as suggested by 
Dalitz et al. [111]. A shell- model analysis of ^"^C seems to support the attractive P-wave AA^ interaction 
[112]. In our QM interaction, this attraction originates from the strong AN-HN^I = 1/2) ^Pi - "^Pi 
couphng due to the FB LS^~^ force [41]. 

Some comments are in order with respect to particle-basis calculations of the Ap and An scatterings. 
From the energy spectrum of the - "^He isodoublet A-hypernuclei, it is inferred that the Ap interaction 
is more attractive than the An interaction. This CSB may have its origin in the different threshold 
energies of the EA^ particle channels and the Coulomb attraction in the E~p channel for the An system. 
The former effect increases the cross sections of the Ap interaction and the latter the An interaction. 
However, the full calculation including the pion- Coulomb correction using the correct threshold energies 
yields very small difference in the S'-wave phase shifts. In the energy region up to pa = 200 McV/c. the 
Ap phase shift is more attractive than the An phase shift only by less than 0.2° in the ^5*0 state, and by 
less than 0.4° in the ^5*1 state. This can also be seen from the results of the Ap and An effective range 
parameters in the particle basis. On the other hand, a rather large effect of CSB is found in the E~p 
channel as discussed in the next subsection. This is because the CSB effect is enhanced by the strong 
AN-T,N(I = 1/2) channel-coupling effect in the ^Si - ^Di state. 
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3.2.4 S p system 

Since the E-p system is expressed as |S^p) = -^/2/3\EN{I = 1/2)) + ^/T/3\EN{I = 3/2)) in the 
isospin basis, it is important to know first the phase-shift behavior of the SiV(J = 1/2) system. The 
SUs decomposition of the SAr(/ = 1/2) state, that is, SAr(/ = 1/2) = [3(11), - (22)] /\/IO for 
and and EN{I = 1/2) = [(11)„ + (03)] /^ for and (see Table 5), is very useful to know 
the QM prediction for the phase-shift behavior in the isospin basis. Since the (ll)s SUs state is 
completely Pauli forbidden, the HN{I = 1/2) ^5*0 phase shift is strongly repulsive due to the exchange 
kinetic-energy kernel [40]. On the other hand, the SA^(/ = 1/2) '^^i phase shift is expected to be 
attractive similarly to the AA^ ^Si phase shift, in so far as the effect of the FSB is not important. 
(Note that AA^ = [-(11)„ + (03)] /\/2 for ^E and states.) Unfortunately the last condition is 
applicable only to the central force, since the one-pion tensor force introduces considerable complexities 
in the SA^(/ = l/2)-AA^ ^5*1 - ^Di CC problem. The strength of the T.N{I = 1/2) central attraction, 
discussed in the preceding subsection, should therefore be examined carefully after this 5'-wave and 
D-wave channel coupling is properly treated. 

Figure 15 displays the ^So phase shift and the -S-matrix of the AN-1:N{I = 1/2) '^^i - CC 
state. The upper figures are the predictions by fss2, while the lower ones by FSS. The ^So phase shift 
shows very strong Pauli repulsion, similarly to the the SA^(/ = 3/2) phase shift (see Fig. 10(a)). 
The SA^(/ = 1/2) ^Si phase shift predicted by fss2 starts from 180° at p^, = 0, decreases moderately 
down to about 160° in the region ps < 100 MeV/c. Then it rapidly decreases to about 80° at around 




200 400 600 800 1000 200 400 600 800 1000 
P£ (MeV/c) Pz (MeV/c) 



Figure 15: The same as Fig. 14 but for the ^5*0 phase shift and for the ^'-matrix of the AA'-EA^(/ = 1/2) 
- ^Di CC state. In the latter system, the channels are specified by 1: AA^ ^5*1, 2: AA^ ^Di, 3: EA^ 
^Si and 4: EA" ^Di. The upper figures display the results by fss2, while the lower by FSS. 
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ps = 120 MeV/c. Beyond this momentum, a moderate decrease follows again. The situation in FSS 
is rather different. The ^5*1 phase shift predicted by FSS starts from 0° at p^. = 0, and shows a clear 
resonance behavior in the region 100 < ps < 120 MeV/c. The peak value of the phase shift is about 
60°. In spite of the very different behavior of the diagonal phase shifts predicted by fss2 and FSS, the 
(S'-matrices are found to be very similar to each other. This can be seen from the Argand diagram 
showing the energy dependence of the T,N{I = 1/2) ^5*1 S-matrix element 6*33 = rj^se^''^'^'^ . This implies 
that the strength of the central attraction in the EA^(/ = 1/2) channel is almost the same as in fss2 
and FSS. 

The very strong reduction of 7733 = |5'33| in Fig. 15 is related to the large enhancement of the 
transmission coefficients, rjis and ?723- from the I1N{I = 1/2) ^Si {i — 3) to the AiV ^5*1 {i = 1) and AiV 
■^Z^i {i = 2) channels around ps = 120 MeV/c. These transmission coefficients are directly connected 
to the T,~p — > An reaction cross sections (which we call cr(C)), and the driving force for this transition 
is the one-pion tensor force. If we calculate the partial-wave contributions to this reaction cross section 
at ps = 160 MeV/c, the sum of the contributions from the transitions ^5"! — > ^Si and ^5"! ^Di 
amounts to about 120 mb both in fss2 and FSS, which is much larger than the contribution (3.5 mb) 
from the transition ^Sq ^Sq- This is in accordance with the analysis of the Ap system in the preceding 
subsection. Namely, there is a large cusp structure in the ^5*1 phase shift at the T,N threshold, while 
a very small cusp in the ^Sq channel. There is, however, some quantitative difference in the detailed 
feature of this tensor couphng between fss2 and FSS. In fss2, 7713 and 7^23 are almost equal to each 
other in the momentum region where the experimental data exist (110 < p^ < 160 MeV/c), whereas 
^13 < V23 holds in FSS. The details of the cross section of the process C indicate that (j{^Si ^Si) ~ 
a{^Si ^D^) in fss2, while a{^Si ^Si) < a{^Si '^D^) in FSS. 

It is important to take into account the pion-Coulomb correction for more detailed evaluation of 
T,~p cross sections. In Ref. [113], we incorporated the Coulomb force of the T,~p channel correctly in 
the particle basis, but the threshold energies of the TPn and An channels were not treated properly. As 
was discussed in Sec. 2.1, we can deal with the empirical threshold energies and the reduced masses in 
the RGM formalism, keeping the correct effect of the Pauli principle. Since the threshold energies and 
reduced masses are calculated from the baryon masses, these constitute a part of the pion-Coulomb 
correction in the YN interaction. Although the pion-Coulomb correction may not be the whole story 
of the CSB, it is certainly a first step to improve the accuracy of the model predictions calculated in 
the isospin basis. We can expect that the small difference of the threshold energies becomes important 
more and more for the low-energy S~p scattering. In particular, the charge-exchange total cross section 
T,~p — > E°7^ (which we call cr(B)) does not satisfy the correct law in the zero-energy limit, if 
the threshold energies of the T,~p and TPn channels are assumed to be equal. We therefore used the 
prescription [114] to multiply the factor (kf/ki), in order to get o'(B) from the cross section calculated 
by ignoring the difference of the threshold energies. Here ki and kf are the relative momentum in 
the initial and final states, respectively. We will see below that this prescription is not accurate and 
overestimates both a(B) and cr(C). 



3.2.5 The inelastic capture ratio at rest in the low-energy S p scattering 



The largest effect of the pion-Coulomb correction appears in the calculation of the E p inelastic capture 
ratio at rest, tr [113]. This observable is defined by [114] 



1 3 

= - rs=o + - rs=i 



with 



^5=0,1 



(T(5=o,l)(B) 



0-(5=0,1)(B) +CT(5=o,l)(C) 



(3.2) 



This quantity is the ratio of the production rates of the E° and A particles when a E particle is trapped 
in an atomic orbit of the hydrogen and interacts with the proton nucleus. For the accurate evaluation 
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Table 7: E~p zero-energy total cross sections kiOs (mb fm~^) for each spin-state (-S" = and 1) and kiOr 
for reactions B: S^j; TPn and C: S^j; An. The divergent factor Cq is taken out in the Coulomb 
case. The inelastic capture ratio tr at rest and in flight are also given. Two versions of our quark 
model, FSS and fss2, are used. 
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40.9 
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of tr, we first determine the effective range parameters of the low-energy ^'-matrix using the multi- 
channel effective range theory [113]. The calculations are performed by using the particle basis with 
and without the Coulomb force. The scattering-length matrices are employed in Table 7 to calculate 
Tr. Compared to the empirical values tr = 0.33 ± 0.05 [115], 0.474 ± 0.016 [116] and 0.465 ± 0.011 
[117], wc find that tr = 0.442 predicted by fss2 is slightly smaller than the recent values between 0.45 
and 0.49. The contribution from each spin state is also listed in Table 7. Wc find rs=Q ~ 0.9, which 
indicates that the cr(C) is very small in comparison with the cr(B) in the spin-singlet state. On the other 
hand, rs=i ~ 0.29 implies that most of the T,N-AN channel coupling takes place in the spin-triplet 
state. Here again the one-pion tensor force is very important in the EA^(/ = l/2)-AiV ^Si - ^Di CC 
problem. We also find that the effect of the Coulomb force plays a minor role for this ratio [113]. 

The E~p inelastic capture ratio in flight rp, predicted by fss2 in the full calculation, is illustrated in 
Fig. 16. Unhke rR, this quantity is defined directly using the total cross sections (7r(B) = (l/4)cro(B) + 
(3/4)(7i(B) and (7t(C) = (l/4)(7o(C) + (3/4)ai(C): 

ariB) 

This quantity is sensitive to the relative magnitudes of the E~p — > E°n and E~p — > An total cross 
sections. In the momentum region of ps > 100 MeV/c, the rp calculated in the particle basis changes 
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Figure 16: E~p inelastic capture ratio in flight 
rp, predicted by fss2, as a function of the incident 
momentum ps- Calculation is made in the parti- 
cle basis with (solid curve) and without (dashed 
curve) the Coulomb force. 
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Figure 17: Calculated low-energy S^p and scattering total cross sections by fss2, compared with 
experimental data: (a) E+p elastic, (b) S~p elastic, (c) T,~p — >• S°n charge-exchange, (d) E~p —>■ An 
reaction cross sections. Predictions in the particle basis without the Coulomb force (dashed curves), 
and those in the isospin basis (dotted curves) are also shown. The experimental data are taken from 
[119] for (a) and (b), and from [118] for (c) and (d). 



little, irrespective of whether the Coulomb force is included or not. The empirical value of rp averaged 
in the interval = HO - 160 McV/c is rp = 0.47 ± 0.03 [118], as shown in Fig. 16 by a cross. The 
prediction of fss2, rp = 0.42, is too small, which is the same feature as observed in FSS {rp = 0.41) 
[113]. The main reason for this disagreement is that our S~p — > An cross sections are too large. 

3.2.6 YN cross sections 

Figure 17 displays the low-energy S^p and S+p cross sections predicted by fss2. The results of three 
different calculations are shown; they are in the particle basis including the Coulomb force (solid curves), 
in the particle basis without the Coulomb force (dashed curves) , and in the isospin basis (dotted curves) . 
The effect of the correct threshold energies and the Coulomb force is summarized as follows. In the 
E+p scattering, the effect of the repulsive Coulomb force reduces the total cross sections by 11 - 6 mb in 
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Figure 18: Calculated E+p and TTp differential cross sections by fss2 (solid curves) and FSS (dashed 
curves), compared with the experimental angular distributions: (a) E+p elastic scattering at = 
170 MeV/c. (b) The same as (a) but at ps = 450 MeV/c. (c) E~p elastic scattering at ps = 160 MeV/c. 
(d) The same as (c) but at ps = 550 MeV/c. (e) E^p — > An difTerential cross sections at ps = 
160 MeV/c. (f) S~p — > E°n differential cross sections at ps = 550 MeV/c. The experimental data are 
taken from [118] for (a), from [120, 121] for (b), from [119] for (c) and (e), and from [122] for (d). 

the momentum range ps = 140 - 180 MeV/c. On the other hand, the attractive Coulomb force in the 
incident E~p channel increases all the cross sections. An important feature of the present calculation 
is the effect of using the correct TPn threshold energy. It certainly increases the E~p TPn charge- 
exchange cross section, but the prescription to multiply the factor (kf/ki) overestimates this effect. The 
real increase is about 1/2-2/3 of this estimation. Furthermore, we find that this change is accompanied 
with the fairly large decrease of the E~p elastic and E~p — > An reaction cross sections. Apparently, this 
effect is due to the conservation of the total flux. The net effect of the Coulomb and threshold energies 
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becomes almost zero for the T,~p elastic scattering. The charge-exchange reaction cross section largely 
increases, and the S~p — > An reaction cross section decreases moderately . The result of fss2 agrees 
with the experimental data reasonably well, although the E~p — > An total reaction cross sections are 
somewhat too large. 

Figure 18 compares the differential cross sections calculated by fss2 and FSS with experiment. For 
E+p and elastic differential cross sections, the recent data taken at KEK [120, 121, 122] are also 
compared. The agreement between theory and experiment is satisfactory. We display in Fig. 19 the 
total cross sections of the NN and YN scatterings in a wide energy range. The solid curves denote the 
fss2 results, while the dashed curves the FSS ones. The "total" cross sections in the charged channels, 
pp, Tj'^p, and arc calculated by integrating the differential cross sections over the angles from 

cos^inin = 0.5 to COS ^max = —0.5. In the NN total cross sections, the effect of the inelastic channel is 
rather small up to Tiab = 400 MeV. For the YN total cross sections, the pion-Coulomb correction is 
important only for the charged channels and the low-energy E~p reaction cross sections. 

3.2.7 G-matrix calculations 

Figure 20 shows various saturation curves for symmetric nuclear matter calculated using the QTQ 
prescription or the continuous prescription for intermediate spectra. The results with the Paris potential 
[84] and the Bonn B potential [123] are also shown for comparison. Since the short-range part of our 
QM interaction is described mainly by the quark-exchange mechanism, its non-local character is entirely 
different from the phenomenological or V-meson exchange pictures in the standard meson-exchange 
models. In spite of this large difference the saturation point of the QM interaction does not deviate 
much from the Coester band, which is similar to other realistic meson-exchange potentials. 

The momentum dependence of the nuclcon. A, and E s.p. potentials Usiq) obtained with the 
continuous prescription is shown in Fig. 21 at three densities p = 0.5 poi 0.7po and po, where po — 
0.17 fm~^ is the normal nucleon density. (These densities correspond to kp = 1.07, 1.2 and 1.35 fm~^, 
respectively.) The results of the Nijmegen soft-core potential NSC89 [124] calculated by Schulze et al. 
[125] are also shown. The corresponding figures of the s.p. potentials calculated by FSS are given in 
Figs. 2 - 5 of Ref. [45]. We see that the nucleon s.p. potential Un{(1) given by fss2 is very similar to 
that of FSS except for the higher momentum region g > 3 fm~^. As already discussed, too attractive 
behavior of FSS in this momentum region is corrected in fss2, owing to the effect of the momentum- 
dependent Bryan-Scott terms involved in the S-meson and V-meson exchange EMEP's. The saturation 
curve in Fig. 20 shows that this improvement of the s.p. potential in the high-momentum region is 
favorable to move the saturation density to the lower side, as long as the calculation is carried out with 
the continuous prescription. On the other hand, the saturation curve with the QTQ prescription alters 
significantly between FSS and fss2. The fss2 prediction with the QTQ prescription is very similar to 
that of the Bonn model-B potential. It is interesting to note that the fss2 result is rather close to Bonn 
model-C for the deuteron properties (see Table 3), but to model-B for the nuclear saturation properties. 
The model-B has a weaker tensor force than the model-C, which is a favorable feature for the nuclear 
saturation properties. 

Figures 21(b) and 21(c) display the momentum dependence of the A and E s.p. potentials in nuclear 
matter, which are obtained from the QM G-matrices of fss2. We find that the Ui^{q) predicted by fss2 
and FSS (shown in Fig. 3 of [45]) are again very similar for < 2 fm~^. The depth of the Up^{q) at g = 
is 48 McV (fss2) and 46 MeV (FSS) in the continuous prescription. This value is slightly more attractive 
than the value expected from the experimental data of A-hypernuclei [126]. On the other hand, the E 
s.p. potential is repulsive in the QM, reflecting the characteristic repulsion in the "^Si + ^Di channel of 
the isospin 7 = 3/2 state (the Pauh repulsion). The depth of the C/s(?) at = is 8 MeV (fss2) and 20 
MeV (FSS). The repulsive feature of the E s.p. potential is supported by Dabrowski's analysis [127] of 
the {K~,T^^) experimental data at BNL [128]. Recently, inclusive {t^~,K^) spectra corresponding to 
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Figure 19: Calculated NN and YN total cross sections by fss2 (solid curves) and FSS (dashed curves), 
compared with the experimental data. The calculations are made in the particle basis. The experimental 
data are taken from [61] for NN, from [108, 109, 110] for Ap, from [118, 120, 121] for E+p, and from 
[117, 119, 122] for E-p scattering. 



the E formation were measured at KEK [129, 130] with better accuracy. These data are satisfactorily 
reproduced using a repulsive E-nucleus potential whose strength is about 30 MeV [131]. 
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Figure 20: Nuclear matter saturation curves ob- 
tained for fss2 and FSS, together with the results 
of the Paris potential [84] and the Bonn model- 
B (Bonn B) potential [85]. The choice of the 
intermediate spectra is specified by "QTQ" and 
"cont." for the QTQ prescription and the contin- 
uous choice, respectively. The result for the Bonn 
B potential in the continuous choice is taken from 
the non-relativistic calculation in Ref . [123] . 
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The partial wave contributions to U\{q = 0) and UY,{q = 0) in symmetric nuclear matter at kp = 
1.35 fm^"*^ are tabulated in Table 8, together with the results of FSS [45] and NSC89 [125]. For the A 
s.p. potential, the characteristic feature of fss2 in comparison with FSS appears in the less attractive 
^5*0 state and in the more attractive ^Si state. The partial wave contributions of fss2 is very similar 
to those of NSC89, except for the ^Si + ^Di contribution. The extra attraction of fss2, compared to 
NSC89, in the A s.p. potential comes mainly from this channel (15 - 16 MeV). This is probably because 
the tensor coupling is stronger in fss2 than in NSC89. A minor excess of the attraction comes from the 
^Pi + ^Pi and ^P2 + ^F2 states (2 - 3 MeV). In Uj:{q = 0), the reduction of repulsion from 20 MeV in FSS 

Table 8: A and S s.p. potentials in nuclear matter with kp = 1.35 fm~^, calculated from fss2 (FSS 
[45]) G-matrices in the continuous prescription for intermediate spectra. Predictions of the Nijmegen 
soft-core potential NSC89 [124] are also shown for comparison [125]. 
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Figure 21: (a) The nucleoli s.p. potential UN{q) in nuclear matter in the continuous choice for inter- 
mediate spectra. Predictions by fss2 for three densities p = 0.5po, 0.7po and po are shown. Here the 
normal density po corresponds to kp — 1-35 fm~^. The dashed curve is the result achieved by Schulze 
et al. [125] with the Nijmegen soft-core NN potential NSC89 [124]. (b) The same as (a) but for the A 
s.p. potential Ui^i^q). (c) The same as (a) but for the E s.p. potential U-E^q). (d) The same as (b) and 
(c), but for the imaginary part of the A and S s.p. potentials Usiq) predicted by fss2. 



to 8 MeV in fss2 is mainly brought about by the 7 MeV reduction of the 1 = 3/2 ^Si + ^Di repulsion 
and by the 4 MeV increase of the 7=1/2 ^5*1 -|- ^Di attraction. The latter feature is again related 
to the strong tensor coupling in fss2. On the other hand, the repulsive contribution of the 7 = 3/2 
^Si + ^Di state in NSC89 is very weak, since this channel has a broad resonance around ps = 500 - 
800 MeV/c. (See Fig. 1 of [99].) It is interesting to note that the attractive contributions to the A and 
E s.p. potentials from the 7 = 1/2 ^51-1-^7^1 state is more than 10 MeV stronger in fss2 than in NSC89. 
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Table 9: Decomposition of the Scheerbaum factors obtained by fss2 at kp — 1-35 fm ^ into various 
contributions. The values in parentheses are given by FSS. The unit is MeV fm^. 
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Although NSC89 is considered to have a strong AN-HN couphng, the AN-IIN coupling of fss2 is even 
stronger. The imaginary parts of the A and E s.p. potentials are shown in Fig. 21(d) for fss2. These 
results are similar to those of FSS (see Fig. 4 of ref. [45]). In particular, Im U^{q = 0) at k^.^ = 1.35 
fm~^ is —14.4 MeV in fss2 and —18.4 MeV in FSS. These results are in accord with the calculations by 
Schulze et al. [125] for NSC89. 

Using the G-matrix solution of fss2, we can calculate the Scheerbaum factor Sb [63], which represents 
the strength of the s.p. spin-orbit potential defined as [46] 

U-(r) = -I ^ . <r . (3.4) 

z r ar 

The value of Sb is momentum-dependent, but here we only discuss (g = 0). The QM description of 
the YN interaction contains the antisymmetric spin-orbit {LS^~^) component, and the large cancellation 
between the LS and LS**^"^ contributions in the AN isospin 7=1/2 channel leads to a small s.p. spin- 
orbit potential for the A-hypernuclei. The FSS gives a ratio Sa/Sn < 1/10 [46]. The cancellation 
of these two contributions is less prominent in fss2, since the S-meson EMEP yields the ordinary LS 
component but no LS^^^ component (see Sec. 2.2). Since the total strength of the LS force is fixed in 
the NN scattering, the FB contribution of the LS force is somewhat reduced. This can easily be seen 
from the simple formula given in Eq. (52) of [46], which shows that in the Born approximation the FB 
LS contribution to the Scheerbaum factor is determined only by a single strength factor asX^mudC%^ 
with X = (h/rriudch). The value of this factor is 29.35 MeV fm^ for fss2, which is just 3/5 of the FSS 
value 48.91 MeV fm^. The density dependence of Sb is found to be rather weak. At the normal density 
po with kp = 1.35 fm~\ we obtain Sn ^ -42.4, Sa = -11.1, Sj^ = -23.3 (MeV hn^) for fss2, resulting 
in Sa/Sn ~ 0.26 and S-z/Sn ~ 0.55. The ratios Sa/Sn and S-^/Sn become slightly smaller for lower 
densities. Table 9 lists the contributions from the LS and LS**-"-* components in the even- and odd-parity 
states as well as / = 1/2 and / = 3/2 channels for kp = 1.35 fm^"^. The parenthesized figures are the 
predictions by FSS. A prominent difference between fss2 and FSS appears only in the contribution 
of Sa- Namely, the 5 MeV reduction of the LS^~^ contribution and the 2 MeV enhancement of the 
LS contribution, which explains the change from Sa = —3.5 in FSS to Sa — —11.1 in fss2. In the 
recent experiment at BNL, very small spin-orbit splitting is reported in the energy spectra of ^Be and 
[132, 133]. A theoretical calculation of these AA^ spin-orbit splittings using OBEP AA^ interactions 
is carried out in Ref. [134]. More detailed analysis of the aaA system for ^Be is made in Sec. 3.5.1, in 
the three-cluster Faddeev formahsm using our QM AA^ interaction [135, 136]. 
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3.3 Interactions for other octet-baryons 
3.3.1 S = —2 systems 
[AA system] 

Although the direct measurement of AA scattering is not possible, the observation of double-A 
hypernuclei gives an important information on the AA interaction. Several candidates of double-A 
hypernuclei are already found [137, 138, 139, 140]. Among them, the sequential decay process of /^jjie 
(NAGARA event), found in the emulsion/scintillating- fiber hybrid experiment at KEK, has given a 
definite value for the binding energy of the two-A hyperons, Saa, without ambiguities arising from the 
possibilities of excited states. The AA interaction energy extracted from Ai?AA(AAHc) = i?AA(AAHc) — 
25A(lHe) is AEaaIaaHc) = 1.01±0.20;2;}? MeV [137]. This value is considerably smaU compared with 
a similar measure of the AA^ interaction strength, A^AArdHe) = SATvdHe) — B\{^Z) — BnC^Rb) — 
^AdHe) — B\{^Z) — 1.73±0.04 MeV, where Ba{^Z) stands for a proper spin and charge average over 
the binding energies of the 0"*" and l"*" states in |H and ^He [141]. 

We display in Figs. 22 and 23 the ""^^o phase shifts of the AA-SA^-SS system with isospin / = 0, 
calculated by FSS and fss2, respectively. The AA-SA^ CC effect gives a step-like resonance (FSS) and 
a cusp (fss2) at the EN threshold region. The couphng effect to the EE channel also gives cusps at the 
EE threshold. The AA^^o phase shift rises up to about 40° for FSS and 18° for fss2. The AA T-matrix 
of FSS yields ABaa ~ 3.7 MeV and that of fss2 ^ 1.4 MeV, as will be shown in the Faddeev calculations 
of the aaHg system. (See Table 17.) Another estimate of AB\\{^fJie), using the CC G-matrices by 
fss2, yields Ai?AA = 1-12 - 1.24 MeV, if the Brueckner rearrangement energy of the a-cluster is properly 
taken into account [142]. These calculations indicate that fss2 is more appropriate than FSS, at least 
for the strength of the AA attraction in j^jjie. 

It is important to understand the origin of the attraction of the ^5*0 A A interaction. In the single- 
channel calculation, it is convenient to display this interaction in terms of the effective local potential 
used in Refs. [97, 35, 98], based on the WKB-RGM approximation. Figure 24 shows the phase-shift 
equivalent local potentials at piab = 200 MeV/c for the ^-S'o states of A^A^, AA^ and AA systems, obtained 
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Figure 22: ^5'o(/ = 0) phase shifts, predicted by 
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Figure 23: The same as Fig. 22, but for fss2. The 
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dicted by the SB effective potential in Eq. (3.12), 
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by solving the transcendental equation of the single-channel RGM kernel of fss2. The single- channel ^5*0 
phase shifts given by these local potentials are depicted in Fig. 25. In the low-energy region, our QM 
gives the attractive A A and AN interactions which are similar to the NN one, in accordance with the 
discussion based on the Hmor-SUs symmetry [98] . The different strengths of the attraction in the NN, 
AN and AA systems originate from the S-meson nonet exchange, i.e., the SU3 relations of the spin- flavor 
factors. On the other hand, the short-range repulsion arising mainly from the color- magnetic term of the 
FB interaction receives a strong effect of the FSB through the mass ratio of the s- and wd-quarks. Since 
the color- magnetic term contains the quark- mass dependence, l/{mimj), the short-range repulsion is 
diminished drastically if the mass ratio A = {iris /mud) deviates largely from unity. Consequently, the 
characteristics of the potential is determined by the subtle balance of the two reductions; one by the 
FSB in the short-range repulsion generated from the color-magnetic term and the other by the SU3 
relations in the medium-range attraction generated from the S-meson nonet exchange. The reduction of 
the attraction from the NN to AN channels suggests that the reduction of the medium-range attraction 
is larger than that of the short-range repulsion in fss2. This feature is also reflected when we extend the 
present framework from the AN to AA systems. The relative strength of the attraction for the diagonal 
^Sq potentials of the NN, AN and AA systems follows the relationship 



U^AA < \Van\ < \Vnn\ 



(3.5) 



This feature for the two-baryon systems including the A-particlc is also seen in NSC97 [143]. 

The full ^Sq potentials are obtained by taking into account the effect of the AiV-SiV channel coupling 
for the AN system and the AA-SA'"-EE coupling for the AA system. First let us discuss the CC effect 
in the ^Sq state without the EMEP. The exchange kinetic-energy, color- Coulombic and momentum- 
dependent Breit retardation terms do not contribute to the AA-EN ^Sq transition potential, because 
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Figure 24: The ^5*0 phase-shift equivalent local 
potentials of single-channel A^A^ (solid curve), AA^ 
(dashed curve) and AA (dotted curve) systems, 
predicted by fss2 at piab = 200 MeV/c. 



Figure 25: A^A^, AA^ and AA ^Sq phase shifts given 
by the phase-shift equivalent local potentials in 
Fig. 24 in the single-channel calculation. 
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the spin-flavor-color factors of these terms are exactly zero, as well as that of the exchange normalization 
kernel. Therefore the contribution from the FB interaction to this coupling potential is solely from the 
color-magnetic term. This coupling potential gives the sharp cusp structure for the AA phase shift 
and the step-like resonance for the EN phase shift, as seen in Fig. 19 of Ref. [11]. If the couphng to 
the EE channel is taken into account, the resonance in the 'EN channel grows up due to the strong 
attraction in the EE channel and moves to the AA channel. The actual resonance behavior is, however, 
sensitive to the FSB. If one neglects the FSB and set A = 1, the CC calculation of FSS without EMEP 
yields a bound state of the binding energy 31 MeV, which corresponds to the iJ-dibaryon predicted by 
Ref. [144]. An introduction of the large FSB reduces the AA-EN CC effect. 

Next, we introduce EMEP's acting between quarks, and investigate how they influence the features 
of the channel coupling. The strange-meson exchange contributes to the AA-SA^ diagram through 
the direct and exchange kernels in the RGM framework, whereas the non-strange-meson exchanges 
can contribute only through the exchange kernel. As the result, the dominant contribution from the 
EMEP's to the AA-EN coupling potential is from the strange-meson exchange. This contribution has 
a sign opposite to that of the color-magnetic term [98, 145] and diminishes the coupling effect by the 
FB interaction discussed above. It turns out that the strange S meson, k, gives the most important 
contribution to the AA-EN ^Sq coupling potential. The effect of the K-meson exchange becomes more 
important as the mass becomes smaller. It should also be stressed that the spin-flavor factors of EMEP's 
are explicitly calculated in our QM baryon-baryon interactions. The contribution from the exchange 
RGM kernel of the EMEP's to the AA-EN ^Sq coupling matrix elements is also appreciable, as largely 
as about one third of the corresponding direct kernel in fss2. Namely, the following three features of 
the EMEP reduce the AA-SA^ CC effect in the present framework; 1) the introduction of the S-meson 
nonet including the k meson, 2) the comparatively small K-meson mass, and 3) the introduction of the 
EMEP's not at the baryon level but at the quark level. 

The above discussion explains why the AA-SA^ CC effect of fss2 is smaller than other several QM's 
[146, 144, 147, 148]. It is related with the more realistic feature including the S-meson nonet (with 
m« = 936 MeV/c^) and the large FSB (with A = 1.5512). After all, the AA interaction of fss2 with 
the full CC effect satisfies the relationship of Eq. (3.5). This feature of the AA interaction leads to 
the conclusion that no bound state is possible in fss2. (For FSS, the relationship Eq. (3.5) is modified 
to \v\\\ ~ \v\n\ < l^iVArl, but it is still insufficient to produce a bound state.) Clearly a realistic QM 
baryon-baryon interaction constructed by supplementing EMEP's to the FB interaction can no longer 
predict a deeply bound i7-dibaryon state. 

In Ref. [149], the AA interaction is examined in the chiral constituent quark model, which is tightly 
constrained by the description of the non-strange sector and the strange meson and baryon spectra. 
They found an attractive AA interaction due to the mechanism that the short-range repulsion from the 
Goldstone-boson and gluon exchanges are considerably reduced in the AA interaction in comparison with 
the A^A^ interaction. The reduction of the short-range repulsion produced from the Goldstone-boson 
exchange might have some relevance in the systems with more strangeness. 

[HiV interaction] 

The SU3 classification in Table 5 is useful to discuss the characteristics of the EN interaction. Let us 
first discuss the phase-shift behavior of the fiavor-symmetric ^Sq state in the single-channel calculation. 
The color-magnetic term is repulsive in the (22) state and attractive in the fiavor-singlet (00) state. 
Since the weights of the (22) and (00) components in the EN {I = 0) state are comparable, the net 
contribution of the color-magnetic term is small in the / = EN interaction due to the cancellation 
between these two components. The EN {I — 0) ^5*0 interaction is, therefore, characterized by the 
EMEP contribution. It turns out that the EMEP contribution produces an attractive EN(I — 0) ^5*0 
interaction. The situation is different in the I — 1 state. In this channel, the leading SU3 component 
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is (ll)s, for which the most compact (Os)® configuration gives a Pauh-forbidden state. This gives a 
repulsive interaction in the SiV(/ = 1) ^Sq channel. As to the flavor-antisymmetric ^Si - ^Di state, the 
EN{I — 0) interaction is found to be weakly attractive, characterizing the property of the pure SU3 
(ll)a state. The EN{I — 1) interaction is repulsive because of the almost-forbidden (Os)^ configuration 
of the (30) state. It is concluded that the SN interaction is attractive in the 1 = channel, whereas it 
is repulsive in the 1 = 1 channel [98]. The attraction of the 5A^(/ = 0) ^5*0 state is strongest, and the 
repulsion of EN{I = 1) ^Sq state is strongest. 

The CC effect on the AA and EE configurations takes palce in the / = ^-S'o state, while that on 
EA in the 1=1 state. In the former case, the attractive feature of the EN interaction does not change 
essentially. On the other hand, the effect of the channel coupling on the SA configuration with 1 = 1 
changes rather drastically the result of the single-channel calculation. Figure 26 displays the phase shifts 
of the ^5*0 and ^5*1 states, calculated for the full 5A^-EA-ES CC system with 1=1. In both cases, the 
phase shifts exhibit a prominent cusp structure at the EA threshold around piab = 600 MeV/c and they 
almost vanish below the EA threshold. This drastic change is due to the strong CC effect caused by 
the large SiV-EA coupling matrix elements. Although the effect of SA^-SA coupling predicted by the 
FB interaction is moderate [10], the introduction of the EMEP's reinforces it since the strange-meson 
exchange can contribute to the SA^-EA coupling through both of the direct and exchange Feynman 
diagrams [150]. 

Figure 27 shows the E~p and E^p total cross sections predicted by FSS (left) and fss2 (right). The 
calculation is made in the particle basis, including the pion-Coulomb corrections. Both of the isospin / = 
and 1 channels contribute to the E~p cross sections, while only the channel with 1 = 1 contributes to 
(or E~n). The E^p (and E~n) total cross sections with the pure 1 = 1 component are predicted to be 
very small below the EA threshold. This behavior of the E~n total cross sections is essentially the same 
as the Nijmegen result [143] . On the other hand, the E~p total cross sections exhibit a typical channel- 
coupling behavior which is similar to that of the E~p total cross sections. These features demonstrate 
that the SA channel coupling is very important for the correct description of scattering observables, 
resulting in the strong isospin dependence of the EN interaction. It is interesting to find that the 
contribution of the strange-meson exchange reduces the AA(7 = 0)-EN{I = 0) channel coupling, 
while it strengthens the EN{I = 1)-EA(7 = 1) coupling in the QM including the FB interaction 
and the EMEP's. This is quite different from the OBEP models, in which both AA-SA" and EN- 
EA coupling potentials are generated only from the strange-meson exchanges. Figure 27 also includes 
the in-medium experimental E~N total cross section around piab = 550 MeV/c [151], where a^-N (in 
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Figure 26: (a) ^5*0 phase shifts in the SA^-EA CC system with 1 = 1, predicted by fss2. (b) ^Si phase 
shifts in the 5A"-EA-EE CC system with 1=1, predicted by fss2. 
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medium) = 30 ± mb is given. Another analysis using the eikonal approximation [152] gives 

a = 20.9 ± 4.5^2 4 '^b, as shown by the black square in Fig. 27. The analysis also estimates the cross 
section ratio a^-p/cTs-n = l-1^0 7-0 4 ^^' ^ ^lore recent experimental analysis [153] for the low-energy 
S~p elastic and S~p — > AA total cross sections in the range of 0.2 GeV/c to 0.8 GeV/c shows that the 
former is less than 24 mb at 90% confidence level and the latter of the order of several mb, respectively. 
These results seem to favor the predictions by fss2. However, we definitely need more experimental 
data with high statistics. A theoretical evaluation of the in-medium cross sections is also necessary. 

The P-wave phase shifts of the AA(/ = 0), EN {I = 0) and 3N{I = 1) scatterings are illustrated 
in Fig. 28 for FSS and in Fig. 29 for fss2. A step-like resonance appears in the ^Pq states of the AA 
channel (FSS) and the EN {I — 0) channel (fss2). The origin of these resonances is the EN (I — 0) 
^Pq resonance produced by the non-central forces, i.e., the tensor and LS forces. This resonance is not 
clearly seen in the E~p total cross sections in Fig. 27, due to the small statistical weight of J = 0. The 
shift of the EN [I = 0) resonance to the AA channel for FSS resembles the situation in the SA^(/ = 1/2) 
^Pi resonance, discussed in Sec. 3.2.3. If the step-like resonance in the T,N{I = 1/2) ^Pi channel moves 
to the AA^ ^Pi channel by the AN-T,N channel couphng, the Ap total cross section is enhanced in 
the cusp region [37, 40]. Similarly, the step- like resonance in the AA "^Pq channel due to the AA- 
EN channel coupling leads to the enhanced AA total cross sections by FSS in the momentum region 
Piab ~ 400 MeV/c. The total cross section curves for the AA scattering clearly show this situation, as 
displayed in Fig. 30. 

[SS interaction] 




Figure 27: E p and S°p total cross sections predicted by FSS (left panels) and fss2 (right panels). The 
calculation is made in the particle basis, including the pion-Coulomb corrections. 
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The properties of the single-channel EE interactions are extensively studied in Ref. [98]. The 
EE(/ = 0) interaction in the ^5*0 state is strongly attractive due to the color-magnetic interaction. 
The central force of the EE(/ = 1) interaction in the ^5*1 state is repulsive due to the Pauli effect. 
Since these 7 = and 1 EE configurations couple to other two-baryon continua, a direct experimental 
confirmation of these single-channel properties is not easy. On the other hand, the 1 — 2 EE channel 
does not couple to any other two-baryon channels. Since the EE(/ = 2) state consists of the pure 
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Figure 28: (left) The P-wave phase shifts of the 
(a) AA(J = 0), (b) EN{I = 0) and (c) EN{I = 1) 
scatterings, predicted by FSS. 
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Figure 30: Total cross sections for the AA scattering, predicted by FSS (a) and by fss2 (b). 



(22) state, the behavior of this interaction is similar to that of the A^A^ ^5*0 state. Figure 9(a) shows 
that this is indeed the case for fss2 and the ^5*0 phase shifts of the A^A^ and SS(/ = 2) channels are 
very similar. Consequently, the and pp total cross sections in Fig. 9(b) are almost equal. 

3.3.2 S = -3 systems 
[HA interaction] 

As discussed in Sec. 3.2.1, the SA-SS interactions are most easily understood on the analogy of 
the AN-TjN interactions. The replacement of A^ with S induces the transition from the (03) to (30) 
SU3 symmetry. This implies that the attractive feature of the spin-triplet deuteron channel is no longer 
taken over to the SA-SE(/ = 1/2) CC system. Figures 31(a) and 31(b) illustrate this situation clearly 
for the ^5*0 and phase shifts predicted by fss2. For the flavor-symmetric ^Sq state, the SU^ contents 
of the SA-SS(/ = 1/2) state are exactly the same as those of the AN-TiN{I = 1/2). Thus the SA and 
SS(/ = 1/2) phase shifts are very similar to the AA^ and UN {I = 1/2) phase shifts, respectively, with a 
slight decrease of the magnitude due to the FSB. On the other hand, the flavor- antisymmetric ^Si state 
of the SA and SE(7 = 1/2) channels yields very similar repulsive phase shifts, related to the strong 
repulsion of the pure (30) UN (I = 3/2) interaction. The fact that the magnitude of the repulsion is 
almost half of the latter is consistent with the weakly attractive nature of the SU^ (ll)a SA^(J = 0) ^5*1 
interaction (see Fig. 10(a)). The EA-EH channel coupling gives a very small cusp at the SE threshold 
both in the ^Sq and ^Si states. It should be stressed that, in the ^Si state, the SA interaction does not 
resemble the AA^ interaction even qualitatively, since the AA^ ^Si state consists of the (ll)a and (03) 
components, while the SA ^Si state consists of (ll)a and (30) components. 

It is reported through Faddeev calculations of light SA hypcrmiclci that =^Hc makes the onset of 
nuclear stability for S hyperons [154]. These authors, however, used an attractive SA effective force 
which simulates the strongly attractive SA interaction predicted by the Nijmegen soft-core NSC97 
model [155, 143]. In view of the relationship between the SA and AA^ channels for ^5*0 and of the 
repulsive SA interaction for ^Si, it is very unlikely that fss2 predicts saHg as a hghtest particle-stable 
S = —3 hypernucleus. 

[HE interaction] 
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As discussed above, the SS interactions with J = 1/2 are repulsive for both ^5*0 and ^5*1 states for the 
kinematical reason. The former is due to the (ll)s Pauh-forbidden state and the latter due to the almost- 
forbidden (30) component. One cannot expect the strong SA-SE coupling by the one-pion exchange 
tensor force unlike the AN-T,N channel couphng, since the SU3 relation of the meson-baryon couphng 
yields Jeett/ Jnntt = "^otm — 1 = for the magnetic-type pseudoscalar coupling with = 2/5 (the 

S'C/e relation). These features are clearly seen in the phase shifts drawn in Figs. 31(a) and 31(b). 

On the other hand, the SS interactions with 1 = 3/2 are attractive for both ^5*0 and states. 
Table 5 shows that the ^-S'o state has the pure (22) symmetry, shared with the NN ^Sq state, and the ^5"! 
state has the (03) symmetry of the NN ^Si state. Of prime importance are the FSB and the fact that 
the SE system is composed of two non-identical particles. For the pion-cxchange contributions, the SUs 
relations of the coupling constants are also important. The SE(/ = 3/2) ^Sq phase shift in Fig. 9(a) 
shows that the effect of the FSB is fairly large for fss2 and the phase-shift rise is about 40 - 50°. On the 
other hand, the similar magnitude of the attraction of the SE(/ = 3/2) '^^i phase shift (Fig. 10(a)) is 
mainly related to the role of the exchange Feynman diagram, which is entirely different from the NN 
system. Unlike the (03) state in A^A^, the (03) state in SE(/ = 3/2) is only moderately attractive, since 
the SE system does not allow the strong one-pion exchange in the exchange Feynman diagram. The 
reduction of the attraction in the direct Feynman diagram is also expected from the SU3 relations for 
the SStt coupling. In summary, the SE interaction predicted by the QM baryon-baryon interaction 
has a large isospin dependence like the EA^ interaction, namely, the repulsive behavior in the 1 — 1/2 
state and the attractive behavior in the 1 — 3/2 state. In particular, fss2 predicts no bound state in 
the S~E~ ^5*0 state. 
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Figure 31: (a) ^S'o phase shifts for the AA^, EA^(/ = 
by fss2. (b) The same as (a) but for the ^5*1 state. 
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3.3.3 S = —4 systems 

The SU3 classification of the SS system shows that the 1 = system is represented by the pure (30) 
symmetry, while the 1=1 system by the pure (22) symmetry. We see in Figs. 10(a) and 9(a) that 
the SS(/ = 0) interaction is repulsive like the EA^(/ = 3/2) one in the ^5*1 state and the SS(/ = 1) 
interaction is attractive like the NN(I = 1) one in the ^-S'o state. In both cases, the effect of the FSB 
is appreciable. The prediction by fss2 gives no bound state in the S°S° ^5*0 state. The SS total cross 
sections are not so large as the NN ones, as seen in Fig. 9(b). 

3.4 Triton and hypertriton Faddeev calculations 
3.4.1 Three-nucleon bound state 

Since our QM EgBg interaction describes the short-range repulsion very differently from the meson- 
exchange potentials, it is interesting to examine the three-nucleon system predicted by fss2 and FSS. 
Here we solve the Faddeev equation for in the framework discussed in Sec. 2.4 directly using the 
QM RGM kernel in the isospin basis [156]. For the NN system, we have no Pauli- forbidden state. 
However, the energy dependence of the two-cluster RGM kernel should be treated properly through the 
self-consistency condition (2.68). Our Faddeev calculation is the full 50-channel calculation, including 
all the angular momentum up to J = 6 for the A'"A^ interaction, and the triton energy converges almost 
completely as seen in Table 10 [157]. The fss2 prediction, —8.52 MeV, seems to be overbound compared 
to the experimental value i?'^^P(^H) = —8.482 MeV. In fact, this is not the case, since the present 
calculation is carried out using the np interaction. The charge dependence of the A^A^ interaction 
indicates that the interaction of the nn system is less attractive than that of the np system. The 
effect of the charge dependence is estimated to be —0.19 MeV for the triton binding energy [85, 158]. 
If we take this into account, our result is still 0.15 MeV less bound. 

For a realistic calculation of the '^H binding energy, it is essential to use such an NN interaction that 
reproduces both the correct D-state probability Pd of the deuteron and the effective range parameters 
of the ^So scattering [159]. Since all the reahstic A'^A'^ interactions reproduce the A^A^ phase shifts more 
or less correctly, the strength of the central attraction is counterbalanced with that of the tensor force. 
Namely, if the interaction has a weaker tensor force, then it should have a stronger central attraction. 
Generally speaking, the effect of the tensor force is reduced in the nuclear many-body systems compared 
to the bare two-nucleon collision. Namely, the A^A^ interaction with a weaker tensor force is favorable 
in order to obtain sufficient binding energies of the nuclear many-body systems. The weak tensor 
force, however, causes various problems such as a too small deuteron quadrupole moment Qa and the 
disagreement of the mixing parameter si of the ^Si + ^Di coupling. 

The correlation between the triton binding energy, Bt = — -E(^H), and the D-state probability of the 
deuteron is plotted in Fig. 32 for fss2 and FSS as well as many other Faddeev calculations using modern 
realistic meson-exchange potentials. In the traditional meson-exchange potentials, the calculated points 
are located on a straight line, which is similar to the Coester line for the saturation point of symmetric 
nuclear matter. For example, the Reid soft-core potential (RSC) [160] gives Pd ~ 6.5 % and predicts 
too small ^H binding energy, Bt = 7.35 MeV. A series of the Bonn potentials reproduce the A^A^ 
phase shifts very accurately, but they tend to have a rather weak tensor force [85] . The model C has the 
strongest tensor force, Pd = 5.61 %, yielding Bf = 7.99 MeV. The value Pd becomes smaller for models 
B and A, and the value of Bt becomes larger accordingly. The following results are given in Ref. [85]: 
model-B (Pd = 5.0 %, Bt = 8.13 MeV), model-A {Pd = 4.4 %, Bt = 8.32 MeV). These results are 
all obtained in the 34 channel calculations (including the A^A^ total angular momentum J < 4), and by 
using the np interaction. In fact, the effects of the charge dependence and the charge asymmetry are 
important for the detailed discussion. The most recent Faddeev calculation employing the CD-Bonn 
potential [47] incorporates these effects, and predicts Bt = 8.014 MeV [162] for Pd = 4.85 %. The 



59 



exp. 



fss2 



Bonn-A 



Bonn-B 



Chiral 



FSS 
Bonn-C 



(Idaho-A) CD-Bonn * 



AV18 



5 6 

Pd (%) 



RSC 



Figure 32: Calculated binding energies Bt as a 
function of the deuteron D-state probability Pp. 
All the calculations are made in the isospin basis, 
using the np interaction, except for the CD-Bonn 
potential. The calculated values are taken from 
[85] (RSC, Bonn-A, B, C), from [163] (AV18 [161], 
CD-Bonn [47]), and from [164] (Chiral). The ex- 
perimental value, Bt = 8.482 MeV, is shown by 
the dashed line. 



present status of the binding energy calculation is summarized as follows: more than 0.5 MeV is 
missing for any realistic two-nucleon interactions [163]. 

On the other hand, our result of the Po value is 5.49 % for fss2 and 5.86 % for FSS. The deuteron 
quadrupole moment Qd is 0.270 fm^ for fss2 and 0.283 fm^ for FSS. Our results of the effective range 
parameters listed in Tables 4 are not as perfect as those of the meson-exchange potentials, but the 
deuteron binding energy and the scattering length as for the ^5*0 state are fitted in determining our 

Table 10: The three-nucleon bound state properties predicted by the Faddeev calculation with fss2 
and FSS. The np interaction is used in the isospin basis. The heading "No. of channels" implies the 
number of two-nucleon channels included; rzmax is the dimension of the diagonalization for the Faddeev 
equation; enn the NN expectation value determined self-consistently; £^(^H) the ground-state energy; 
and \/Jr^)3}i (-\/('"^)3He) the charge rms radius for (^He), including the proton and neutron size 
corrections through Eq. (3.6). The Coulomb force and the relativistic corrections are neglected. 



model 


No. of 


^max 




E(3H) 




V'('^^)3He 




channels 




(MeV) 


(MeV) 


(fm) 


(fm) 




2(3) 


2,100 


2.361 


-7.807 


1.80 


1.96 




5 (SD) 


5,250 


4.341 


-8.189 


1.75 


1.92 


fss2 


10 (J < 1) 


10,500 


4.249 


-8.017 


1.76 


1.94 




18 (J < 2) 


18,900 


4.460 


-8.439 


1.72 


1.90 




34 (J < 4) 


35,700 


4.488 


-8.514 


1.72 


1.90 




50 (J < 6) 


112,500 


4.492 


-8.519 


1.72 


1.90 




2(^) 


2,100 


2.038 


-7.675 


1.83 


1.99 




5 (SD) 


5,250 


3.999 


-8.034 


1.78 


1.95 


FSS 


10 (J < 1) 


10,500 


3.934 


-7.909 


1.78 


1.97 




18 (J < 2) 


18,900 


4.160 


-8.342 


1.74 


1.93 




34 (J < 4) 


35,700 


4.175 


-8.390 


1.74 


1.92 




50 (J < 6) 


112,500 


4.177 


-8.394 


1.74 


1.92 
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QM parameters in the isospin basis. Under these circumstances, the results of fss2 and FSS in Fig. 32 
seem to indicate that they form another - Pjj hne deviating from the traditional meson-exchange 
potentials. After the correction by the charge dependence of the NN interaction, our binding 
energies are smaller than the experimental value by merely 0.15 MeV for fss and 0.28 MeV for FSS. 
If we attribute these differences to the effect of the three-nucleon force, it is by far smaller than the 
generally accepted values, 0.5 - 1 MeV [163], obtained with the meson-theoretical A^A^ interactions. A 
remarkable point of our results is that we can reproduce enough binding energy of the triton without 
reducing the deuteron D-state probability. 

Table 10 lists the charge rms radii of and ^He, predicted by fss2 and FSS. The finite size corrections 
of the nucleons are made through 

(r2)3H = [i?c(^H)]^ + (0.8502)2 _ g ^ (0.3563)^ , 

(r2)3He = [i?c('He)] ' + (0.8502)^ - i x (0.3563)^ , (3.6) 

where Rc^ stands for the squared charge radius for the point nucleons which is obtained from the 
Faddeev calculations. Since our three-nucleon bound state wave functions are given in the momentum 

representation, we first calculate the charge form factors Fc{Q^), according to the formulation given in 
Ref. [165]. The values Rc^ are then extracted from the power scries expansion of Fc{Q'^) with respect 
to Q^. In the present calculation, the Coulomb force and the relativistic correction terms [166] of the 
charge current operator are entirely neglected. The experimental values are difficult to determine, as 
discussed in Ref. [166]. Here we compare our results with two empirical values 

/^^^_ r 1.70 ± 0.05 fm [167] rr^^ — _ / 1.87 ± 0.05 fm [167] , . 

)3H - I ^ g^^QQ5 [168] ' ^('^ )3He- I ^ g3^QQ3^^ pgj . (3.7) 

The agreement with the experiment is satisfactory both for fss2 and FSS. 

For the systems composed of three identical particles, the self-consistent energy of the two-cluster 
RGM kernel, Snn, has a clear physical meaning related to the decomposition of the total triton energy 
E into the kinetic-energy and potential-energy contributions. A simple manipulation of the Faddeev 
equations and the self-consistency condition yields 

{Ho) = 2{?>s„-E) , {V) = 2>{E-2e^) , (3.8) 

where the expectation value is calculated with respect to the three-nucleon total wave function. From the 
first equation, the Enn is obtained from the kinetic-energy contribution through enn = -E/3 + {Ho)/6. 
Table 11 hsts this decomposition, together with the results of CD-Bonn [47] and AV18 [161] potentials 
[163]. We see that our QM results of Sjvat by FSS and fss2 are just in between these potentials, which 



Table 11: Decomposition of the total triton energy E into the kinetic-energy and potential-energy 
contributions: E = (Hq) + (V). The unit is in MeV. In the present framework, this is given by the 
expectation value e^n of the two-cluster Hamiltonian with respect to the Faddeev solution, which is 
determined self-consistently. The results of CD-Bonn [47] and AV18 [161] are taken from Ref. [163]. 



model 


^NN 


E 




(V) 


fss2 


4.492 


-8.519 


43.99 


-52.51 


FSS 


4.177 


-8.394 


41.85 


-50.25 


CD-Bonn 


3.566 


-8.012 


37.42 


-45.43 


AV18 


5.247 


-7.623 


46.73 


-54.35 
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have very different strengths of the tensor force. The energy-dependence of the A^A^ RGM kernel with 
the exphcit structure +eK has some effects favorable to give a larger binding energy of the triton, since 
the value of the exchange normalization kernel K is —1/9 for the most compact (Os) wave function 
between the two nucleon-clusters. The RGM kernel becomes more attractive for a larger value ol Snn- 
Table 11 indicates a larger (y) for a larger Snn, whereas the CD-Bonn and AV18 potentials do not 
have such an energy dependence. 

The Faddeev calculations for using the QM NN potentials have also been carried out by Takeuchi, 
Cheon and Redish [169] and by the Salamanca- Jiihch group [170]. In the former calculation, the model 
QCM-A gives the NN phase shifts with almost the same accuracy as our model FSS, and predicts 
Pd = 5.58 % and Bt = 8.01 - 8.02 MeV in the 5-channel calculation. These are very similar to our 
FSS results. On the other hand, the Salamanca- Jiilich group predicts Bf = 7.72 MeV, in spite of the 
small D-state probability = 4.85 %. It is not clear to us how they treated the energy dependence 
of the RGM kernel when the separable expansion is introduced for solving the Faddeev equations. We 
think that the fit of the NN phase shifts for higher partial waves, especially for the P waves, has to be 
improved in order to extend their calculation to more than 5 channels. 

3.4.2 The hypertriton 

Next, we apply our QM A^A^ and yA^ interactions to the hypertriton (^H) which has the small separation 
energy of the A-particle, B^^"^ = 130 ± 50 keV [107]. Since the A-particle is far apart from the two- 
nucleon subsystem, the on-shell properties of the AA^ and J^N interactions are expected to be well 
reflected in this system. In particular, this system is very useful to learn the relative strength of ^5*0 
and ^5"! attractions of the AA^ interaction, since the ^5*0 component plays a more important role than 
^Si in this system and the available low-energy Ap total cross section data cannot discriminate many 
possible combinations of the ^Sq and ^Si interactions. In fact, Refs. [171, 172, 173] showed that most 

Table 12: Results of the hypertriton Faddeev calculations by fss2 and FSS. The heading E is the 
hypertriton energy measured from the ANN threshold, B\ the A separation energy, enn {san) the 
A^A" (A A") expectation value determined self-consistently, and Ps the T,NN probability in percent. 
The calculated deuteron binding energy is ea = 2.2247 McV for fss2 and 2.2561 MeV for FSS {sf"^ = 
2.224:64:4: ± 0.000046 MeV [86]). The norm of admixed redundant components is less than 10~^. 



model 


No. of 


E 


Ba 










channels 


(MeV) 


(keV) 


(MeV) 


(MeV) 


(%) 




6(3) 


-2.362 


137 


-1.815 


5.548 


0.450 




15 (SD) 


-2.423 


198 


-1.762 


5.729 


0.652 


fss2 


30 (J < 1) 


-2.403 


178 


-1.786 


5.664 


0.615 




54 (J < 2) 


-2.498 


273 


-1.673 


5.974 


0.777 




102 (J < 4) 


-2.513 


288 


-1.658 


6.022 


0.804 




150 (J < 6) 


-2.514 


289 


-1.657 


6.024 


0.805 




6(5) 


-2.910 


653 


-1.309 


3.984 


1.022 




15 (SD) 


-2.967 


710 


-1.433 


6.171 


1.200 


FSS 


30 (J < 1) 


-2.947 


691 


-1.427 


6.143 


1.191 




54 (J < 2) 


-3.121 


865 


-1.323 


6.467 


1.348 




102 (J < 4) 


-3.134 


877 


-1.317 


6.488 


1.360 




150 (J < 6) 


-3.134 


878 


-1.317 


6.488 


1.361 
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Table 13: Correlation between ^5*0 and ^Si effective range parameters of various AN interactions and 
the A separation energy B\ of the hypertriton. The value for NSC89 is taken from Ref. [173]. The 
experimental B\ value is 130 ± 50 keV [126]. 



model 


as (fm) 


Vs (fm) 


at (fm) 


n (fm) 


Ba (keV) 


FSS [40] 


-5.41 


2.26 


-1.02 


4.20 


878 


fss2 [42, 43] 


-2.59 


2.83 


-1.60 


3.01 


289 


NSC89 [124] 


-2.59 


2.90 


-1.38 


3.17 


143 



meson-theoretical interactions fail to bind the hypertriton except for the Nijmegen soft-core potentials 
NSC89 [124], NSC97f and NSC97e [155]. It is also pointed out in Refs. [171] and [173] that a small 
admixture of the TjNN components less than 1% is very important for this binding. We therefore carry 
out the ANN-TiNN CC Faddeev calculation, by properly taking into account the existence of the SU3 
Pauli-forbidden state (ll)s at the quark level. The necessary formulation is given in Sec. 2.4. 

Table 12 shows the results of the Faddeev calculations using fss2 and FSS. In the 15-channel cal- 
culation including the S and D waves of the NN and YN interactions, we have already obtained 
B\ = —Ed — £"(^11) fti 200 keV for fss2. The convergence in the partial-wave expansion is very rapid, 
and the total angular-momentum J < 4 of the baryon pairs is enough to obtain 1 keV accuracy. 
With 150-channel ANN and IINN configurations we obtain = 289 keV with the J^NN component 

= 0.80% for the fss2 prediction, and Ba = 878 keV with = 1.36% for FSS. 

Table 13 shows the correlation between the A separation energy Ba and the ^5*0 and ^Si AN effective 
range parameters of FSS, fss2 and NSC89. Although all of these AA^ interactions reproduce within the 
experimental error bars the low-energy AA^ total cross section data below about pa = 300 MeV/c, our 
QM interactions seem to be slightly more attractive than the NSC89 potential [124]. The model FSS 
gives a large over binding since the ^Sq AN interaction is strongly attractive. The phase-shift difference 
of the ^So and ^Si states at about pa = 200 MeV/c is A5 = 6^80) - 5{^Si) ^ 30° for FSS, while 
A5 ~ 10° for fss2. Since the fss2 result is still slightly ovcrbound, this difference should be made 
smaller in order to reproduce the experimental value. From the two results given by fss2 and FSS, we 
extrapolate the desired difference to be 0° - 2° in our QM. The same conclusion is also obtained in 
Ref. [174], using simulated interactions of the Nijmegen models. 

In order to make sure that this extrapolation gives a good estimation, we modify the K-meson mass of 
the model fss2 from the original value, = 936 MeV (see Table 2), to 1,000 MeV, and repeat the whole 

Table 14: Decomposition of the A^A^ and YN expectation values {enn and Eyn), the deuteron energy 
{—ed) and the total three-body energy E to the kinetic-energy and potential-energy contributions. The 
unit is in McV. The results for NSC89 arc taken from Ref. [171]. 



model 


hjsfN 


+ vnn = 


^NN 


hd + Vd — —Sd (deuteron) 


FSS 


19.986 - 


-21.303 = 


-- -1.317 


16.982 - 19.238 = -2.256 


fss2 


19.371 - 


- 21.029 = 


= -1.657 


17.495 - 19.720 = -2.225 


NSC89 


20.48 - 


- 22.25 = 


-1.77 


19.304- 21.528 = -2.224 


model 


hvN 


+ Vyn = 




(H) + {V) = E 


FSS 


10.036 


- 4.602 = 


= 5.435 


27.372 - 30.506 = -3.134 
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calculation. This modification makes the ^5*0 AA^ interaction less attractive and the more attractive. 
We obtain i?A = 145 keV with = 0.53 %, which is very close to the NSC89 prediction B\ = 143 
keV with Ps = 0.5% [171, 173]. The effective range parameters of this modified fss2 interaction are 
Qs — —2.15 fm, Vs — 3.05 fm, and at — —1.80 fm, = 2.87 fm. The phase-shift difference is only 1.3° 
and the increase of the total cross section of the AA^ scattering is at most 10 mb at p\ = 100 MeV/c 
(from 286 mb to 296 mb), which is still within the experimental error bars. It should be kept in mind 
that the effective range parameters or the 5- wave phase-shift values determined in this way is very much 
model dependent, since the B\ value is not determined by these quantities alone. It depends on how 
higher partial waves contribute and also on the details of the AN-T,N couphng of a particular model. 
The reliability of the extrapolation shown here is based on the fact that the models fss2 and FSS have 
a common framework for the quark sector and EMEP's. 

Table 12 also shows that the expectation value of the NN Hamiltonian, Enn, determined self- 
consistently is close to the deuteron energy —6^, especially in fss2. This feature is even more conspicuous 
if these energies are decomposed to the kinetic-energy and potential-energy contributions. Table 14 
shows this decomposition with respect to fss2, FSS and NSC89. (For this comparison, we use the 
definition of the kinetic-energy part of the deuteron, given by hd = {Xd\hNN\Xd) / {Xd\Xd) , where Xd is 
the RGM relative wave function between the neutron and the proton.) In fss2, the kinetic energy of the 
A"A^ subsystem is 1.88 MeV larger than that of the deuteron, which suggests that the A^A^ subsystem 
shrinks by the effect of the outer A-particle, in comparison with the deuteron in the free space. In 
NSC89, this difference is smaller, i.e., 1.18 MeV. These results are consistent with the fact that the 
hypertriton in NSC89 is more loosely bound {B\ = 143 keV) [173] than in fss2 (289 keV), and the 
A-particle is very far apart from the NN subsystem. Table 14 lists the kinetic-energy and potential- 
energy decompositions for the AN-T,N averaged YN expectation value Eyn and the total energy E. 
The kinetic energies of Eyn are much smaller than those of Satat, which indicates that the relative 
wave function between the hyperon and the nucleon is widely spread in the configuration space. The 
comparison of the total-energy decomposition shows that the wave functions of fss2 and NSC89 may 
be very similar. A clear difference between them, however, appears in the roles of the higher partial 
waves. The energy increase due to the partial waves higher than the S and D waves is 91 keV in fss2 
and 168 keV in FSS, respectively. On the other hand, the results in Refs. [171] and [173] indicate that 
this is only 20 - 30 keV in the case of NSC89. This difference originates from both of the NN and YN 
interactions. Since the characteristics of the meson-theoretical YN interactions in higher partial waves 
are a priori unknown, more detailed analysis of the fss2 results sheds light on the adequacy of the QM 
baryon-baryon interactions. 

3.5 Application to ^Be and ^^He systems 
3.5.1 The cxcxA system for 

As a typical example of three-cluster systems composed of two identical clusters, we apply the present 
formalism to the aaA Faddeev calculation for ^Be, using the aa RGM kernel and the Aa folding 
potential generated from a simple AA^ effective interaction [135, 136]. For the aa RGM kernel, we use 
the three-range Minnesota force [69] with the Majorana exchange mixture u — 0.94687, and the h.o. 
size parameter ly = 0.257 fm~^. The effective AA" central interaction, denoted SB in Table 15, is a 
Minnesota-type potential 



u 2 — u 

- H Pr 

2 2 



(3.9) 



where v{}E) and v{^E) are simple two-range Gaussian potentials generated from the ^5*0 and 
phase shifts predicted by the QM A A" interaction fss2. The inversion method based on supersymmetric 
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Figure 33: AN-T,N ^Sq (a) and ^Si (b) phase shifts for the isospin 1 — 1/2 channel, calculated with 
fss2 (solid and dashed curves) and the single-channel SB potentials (/ = 1) in Eq. (3.10) (circles). 



quantum mechanics is used to derive phase-shift equivalent local potentials [175]. These potentials are 
then fitted by two-range Gaussian functions: 

v{^So) = -128.0 exp(-0.8908 r^) + 1015 exp(-5.383 r^) (McV) , 

v{^Si) = -56.31 /exp(-0.7517r2) + 1072 exp(-13.74r2) (MeV) , (3.10) 

where / is an adjustable parameter and r is the relative distance between A and N measured in fm. 
Figure 33 shows that these potentials fit the ^5*0 and ^Si AN phase shifts obtained by the full AN-T^N 
CCRGM calculations of fss2. In the ^Si state, the phase shifts are fitted only in the low- momentum 
region with piab < 600 MeV/c, since the cusp structure is never reproduced in the single-channel 
calculation. Since any central and single-channel effective AA^ force leads to the well-known ovcrbinding 
problem of ^He by about 2 McV (in the present case, it is 1.63 McV) [111], the attractive part of the 
^5*1 AA^ potential is modified to reproduce the correct binding energy, £^®^P(^He) = —3.12 ± 0.02 MeV 
with / = 0.8923. This overbinding problem is mainly attributed to the Brueckner rearrangement effect 
of the a-cluster, originating from the starting energy dependence of the bare two-nucleon interaction 
due to the addition of an extra A-particle [176]. The odd-state AA^ interaction is introduced in Eq. 
(3.9) with the Majorana exchange parameter u ^ 1. All partial waves up to AMax = ^iMax = 6 for the 
aa and Aa pairs are included. The direct and exchange Coulomb kernel between the two a-clusters is 
introduced at the nucleon level with the cut-off radius, Rc — 10 fm. In the first calculation using only 
the central force, the SB potential with the pure Serber character (it = 1 in Eq. (3.9)) can reproduce 
the energies of the ground and excited states within 100 - 200 keV accuracy [135]. 

As already discussed in Sec. 3.2.7, one of the interesting problems in ^Be is the very small spin-orbit 
(is) splitting, AEes = E^{3/2+) - E^{5/2+), of the 5/2+ and 3/2+ excited states. The recent Hyperball 
7-ray spectroscopy experiment predicts AE'^^'' = 43±5 keV [132, 133]. In the non-relativistic models of 
the YN interaction, this phenomenon is usually described by the strong cancellation of the ordinary LS 
component and the antisymmetric LS component (LS**^"^ force), the latter of which is a characteristic 
feature of baryon-baryon interactions between non-identical baryons. For example, our FSS interaction 
yields such a strong LS^~^ component [29] that is about one half of the ordinary LS component with 
the opposite sign. We performed the G-matrix calculation in symmetric nuclear matter, using this QM 
interaction [45] , and calculated the Scheerbaum factor Sb [46] . The ratio of Sb to the nucleon strength 
S]sf ~ —40 MeVfm^ is S^/Sj^f ~ 1/5 and Sy:/Sn ~ 1/2 in the Born approximation, while it changes to 
Sa/Sn ~ 1/12 in the G-matrix calculation of the model FSS at the normal density po = 0.17 fm~^. 
The significant reduction of Sa in the latter is due to the enhancement of the LS^~^ component in the 
diagonal AA^ channel, owing to the P-wave AN-TjN coupling. 
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Table 15: The ground-state energy £'gr(l/2+), the 5/2+ and 3/2+ excitation energies, £'x(5/2+) and 
£'x(3/2+), and the spin-orbit splitting, AE'^^ = £'x(3/2+) — £'x(5/2+), calculated by solving the Faddeev 
equations for the aaA system in the jj coupling scheme. The exchange mixture parameter of the SB 
AA^ force is assumed to be u = 0.82. The Aa spin-orbit force is generated from the Born kernel of the 
FSS and fss2 AA^ LS interactions. For the fss2 LS interaction, the LS component from EMEP's is also 
included. The AA^ central potentials v^]^ used are Eqs. (3.9) and (3.10) (SB), the simulated Nijmegen 
models (NS, ND, NF) and Jiihch model-A and B potentials (JA, JB) [177]. 
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In the aaA three-cluster model, the origin of the is splitting of the 5/2+ and 3/2+ excited states 
is only from the Aa is interaction. We therefore carry out the Faddeev calculation in the jj-coupling 
scheme, by using the central plus spin-orbit Aa interactions [136]. The Aa spin-orbit interaction is 
generated from the Born kernel of the AA^ LS QM interaction Eq. (2.31) by the a-cluster folding. The 
three possible LS types, LS, LS^~\ and LS^~^a, in Eq. (2.32) are all included. The derivation of the 
Aa Born kernel from the AA^ Born kernel can be carried out analytically, yielding the final form 

VkaiQf.^i) = E [^r V^''{qf,q,)+X^ V^'%qf,q,)] in • 5a . (3.11) 
r 

Here X-^'^ and V^'^'^''^(q/, qi) are the spin-flavor factors and the spatial integrals, respectively, for the QM 
interaction types, T = and S {S'), in the AN interaction. Since the QM AN interaction contains 

the strangeness exchange term, both the direct {d) and exchange (e) Feynman diagrams between A and 
the nucleon in the a-cluster contribute, and they are extremely non-local. The LS components from 
EMEP's are also included for the model fss2. 

Table 15 shows the results of Faddeev calculations in the jj-coupling scheme, predicted by our QM 
AN LS interaction plus the SB and other various AN central potentials [177]. The ground-state energy 
does not change much from the LS'-coupling calculation, indicating the dominant s-wave coupling of the 
A-hyperon. We can reproduce the ground-state energy using the SB central force with the Majorana 
exchange mixture u = 0.82. The final values for the is splitting of the 5/2+ - 3/2+ excited states 
are AE^g = 137 keV (FSS) and 198 kcV (fss2). Compared to the experimental value of 43 ± 5 keV, 
these predictions are three to five times too large. For the G-matrix simulated NSC97f LS potential 
in Ref. [134], we obtain 209 keV for the same SB force. The deviation from 0.16 MeV in Ref. [134] is 
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due to the model dependence of the aa and Ka interactions. Our FSS prediction for AE^s is less than 
2/3 of the NSC97f prediction, while fss2 gives almost the same result as NSC97f. If we switch off the 
EMEP contribution in the fss2 calculation, we find AE^g = 86 keV. This results from the dominant LS 
component (and the very small LyS^"^ component) generated from the EMEP of fss2. 

We see that the direct use of the QM Born kernel for the AA^ LS component is not good enough to 
reproduce the small experimental value 43 ± 5 kcV for the is splitting. In order to clarify the origin of 
the problem, we have carried out an analysis of the Aa s.p. spin-orbit potential, using the Scheerbaum 
factor, [136] . Table 16 lists Sa in symmetric nuclear matter, obtained by the G-matrix calculations 
in the continuous prescription. The Fermi momentum, kp — 1-07 fm~^, corresponding to the half of 
the normal density po = 0.17 fm~^, is assumed. The table also shows the decompositions into various 
contributions and the results obtained by turning off the AA^-SA^ coupling through the LS'^^^ and 
LS^~^a forces. For FSS, we find a large reduction of Sa from the Born value —7.8 MeV fm^ [46], 
especially when the P-wave AN-T,N couphng is properly taken into account. When all the AN-T,N 
couplings including those by the pion tensor force are turned off, the LS^~^ contribution is just a half of 
the LS contribution with the opposite sign (in the dominant odd partial wave), which is the same result 
as in the Born approximation. The + ^-^2 AA^-SA^ coupling slightly enhances the attractive LS 
contribution, while the ^Pi + ^Fi AN-HN coupling largely enhances the repulsive LS^~^ contribution. 
If we use this reduction of the Sa factor from —7.8 MeV fm^ to —1.9 MeV fm^ in the reahstic G- 
matrix calculation, the AE^s value is reduced from —137 MeV to an almost correct value —33 keV. 
However, such a reduction of the Scheerbaum factor due to the AN-HN coupling is supposed to be 
hindered in the Aa system in the lowest-order approximation from the isospin consideration. On the 
other hand, the situation of fss2 in Table 16 is rather different, although the cancellation mechanism 
between the LS and LS^~^ components and the reduction effect of -S'a factor in the full calculation 
are equally observed. When all the AN-T,N couphngs are neglected, the ratio of the LS'^~^ and LS 
contributions in the quark sector is still one half. Since the EMEP contribution is mainly for the LS 
type, it amounts to about —6 MeV fm^, which is very large and remains with the same magnitude 
even after the P-wave AN-HN coupling is included. Furthermore, the increase of the LS^~^ component 
is rather moderate in comparison with the FSS case. This is because the model fss2 contains an 
appreciable EMEP contribution (about 40%) which has very little LS^~^ contribution. As a result, 
the total S\ value in fss2 G-matrix calculation is 3 - 6 times larger than the FSS value, depending 
on the Fermi momentum kp = 1.35 - 1.07 fm~ . Such an appreciable EMEP contribution to the LS 
component of the YN interaction is not favorable to reproduce the negligibly small is splitting of ^Be. 

Table 16: The Scheerbaum factor Sa in symmetric nuclear matter with kp = 1-07 fm~^, predicted by 
G-matrix calculations of FSS and fss2 in the continuous prescription for intermediate spectra. Decom- 
positions into various contributions are shown, together with the cases when the AN-J^N coupling by 
the LS'^'^ and LS'^~^a forces is switched off (couphng off). The unit is in MeV fm^. 
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In conclusion, the failure of the present Faddeev calculation is probably because the P-wave AN-YjN 
coupling is not properly taken into account. The QM baryon-baryon interaction with a large spin-orbit 
contribution from the S-meson exchange potentials is not appropriate to reproduce the very small is 
splitting observed in ^Be. A way of incorporating the P-wave AN-T,N couphng in the cluster model 
calculations is a future problem. 

3.5.2 The AAa system for j^^He 

As another application of the Aa T-matrix in the aaA Faddeev calculation, we discuss the ground-state 
energy of fj^He [178]. The full CC T-matrices of fss2 and FSS with strangeness S = —2 and isospin 
/ = are employed for the AA RGM T-matrix (T-matrix). Table 17 shows the two-A separation energy 
APaa (defined by APaa = Baa{a\^^) ~ 2PA(AHe)), predicted by various combinations of the AN 
and AA interactions. The results of a simple three-range Gaussian potential, i'AA(Hiyama), used in 
Ref. [177] are also shown. This AA potential and our Faddeev calculation using the FSS T-matrix yield 
very similar results with large APaa values of about 3.6 McV, in accordance with the fact that the AA 
phases shifts predicted by these interactions increase up to about 40° (see Fig. 22). The improved QM 
fss2 yields APaa — 1-41 MeV. (For the AA single-channel T-matrix, this value is reduced to 1.14 MeV.) 
In Table 17, results are also shown for vaa(SB) which is generated from the fss2 ^5*0 AA phase shift 
(see Fig. 23) in the full-channel calculation, using the supersymmetric inversion method [175]. This 
single-channel effective potential is given by 

^;aa(SB) = -103.9 exp(-1.176 r^) + 658.2 exp(-5.936 r^) (MeV) . (3.12) 

We think that the 0.5 MeV difference between our fss2 result and the i'aa(SB) result is probably because 
we neglected the CC effects of the AAa channel to SiVa and EEa channels. We should also keep in 
mind that in all of these three-cluster calculations the Brueckner rearrangement effect of the a-cluster is 
very important, producing about —1 MeV (repulsive) contribution [142]. It is also reported in Ref. [179] 
that the quark Pauli effect between the a-cluster and the A-hyperon yields a non-negligible repulsive 
contribution of 0.1 - 0.2 MeV for the two-A separation energy of ^^He, even when a rather compact (3g) 
size of 6 = 0.6 fm is assumed as in our QM interactions. Taking all of these effects into consideration, 
we can conclude that the present results by fss2 are in good agreement with the recent experimental 
value, AB^^l = 1.01 ± 0.20 MeV, deduced from the NAGARA event [137]. 

Table 17: Comparison of APaa values in MeV, predicted by various AA interactions and AA^ potentials. 
The AA potential f AA(Hiyama) is the three-range Gaussian potential used in Ref. [177], and f aa(SB) the 
two-range Gaussian potential given in Eq. (3.12). FSS and fss2 use the AA RGM T-matrix (T-matrix) 
with £aa, which is the AA expectation value determined self-consistently. The AA^ central potentials 
^ATv used are Eqs. (3.9) and (3.10) (SB), the simulated Nijmegen models (NS, ND, NF) and Jiilich 
model-A and B potentials (JA, JB) [177]. The experimental value is AP^\p = 1.01 ± 0.20 MeV [137]. 
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4 Summary 



Since the advent of the Yukawa's meson theory, an enormous amount of efforts have been devoted 
to understand the fundamental micleon-nucleon (NN) interaction and related hadronic interactions. 
The present view of these interactions is a non-perturbative realization of inter-cluster interactions, 
governed by the fundamental theory of the strong interaction, quantum chromo dynamics (QCD), in 
which the gluons are the field quanta exchanged between quarks. On this basis, the meson "theory" 
can be understood as an effective description of the quark-gluon dynamics in the low-energy regime. 
The short-range part of the A^A^ and hyperon-nucleon (YN) interactions are still veiled with unsolved 
mechanism of quark confinement and multi-gluon effects. 

In this review article, we applied a constituent quark model to the study of the baryon-baryon 
interactions, in which some of the essential features of QCD are explicitly taken into account in the 
non-relativistic framework. For example, the color degree of freedom of quarks is explicitly included, and 
the full antisymmetrization of quarks is carried out in the resonating-group method (RGM) formalism. 
The gluon exchange effect is represented in the form of the quark-quark interaction, for which a color 
analogue of the Fermi-Breit (FB) interaction is used with an adjustable parameter of the quark-gluon 
coupling constant as- The confinement potential is a phenomenological r^-type potential, which has no 
contributions to the baryon-baryon interactions in the present framework. Since the meson-exchange 
effects are the non-perturbative aspect of QCD, these are described by the effective-meson exchange 
potentials (EMEP's) acting between quarks. 

Our purpose is to construct a realistic model of the A^A^ and YN interactions, which describes 
not only the baryon-baryon scattering quantitatively in the wide energy region, but also reproduces 
rich phenomena observed in few-baryon systems and various types of finite nuclei and nuclear matter. 
The present framework incorporates both the quark and mesonic degrees of freedom into the model 
explicitly. This framework is very versatile, since it is based on the natural picture that the quarks 
and gluons are the most economical ingredients in the short-range region, while the meson- exchange 
processes dominate in the medium- and long-range regions of the interaction. The short-range repulsion 
is described by the color-magnetic term of the FB interaction and the Pauli effect at the quark level. In 
this article, we have focused on our recent versions called FSS [39, 40, 41] and fss2 [42, 43]. In FSS only 
the scalar- meson and pseudoscalar- meson exchanges are introduced as the EMEP's, while in fss2 the 
vector-meson exchanges are introduced as well. All possible standard terms used in the non-relativistic 
one-boson exchange potentials (OBEP) can in principle be included in the final version. As a first step 
to study the charge dependence of the baryon-baryon interactions, the pion-Coulomb correction is also 
taken into account in the particle basis. After these improvements mainly related to the EMEP's, the 
most recent model fss2 has achieved an accurate description of the A^A^ and YN interactions. 

An advantage of using the quark model for the study of the baryon-baryon interactions is that the 
fiavor 5*^73 symmetry (and also the spin-fiavor 5*^76 symmetry) is automatically incorporated in the 
framework, so that all the interactions between the octet-baryons (Sg) are treated in a unified manner 
in close connection with the well-known A^A^ interaction. We extended the {3q)-{3q) RGM study of 
the the A^A^ and YN interactions to the strangeness S = —2, — 3 and —4 sectors with no ad hoc 
parameters, and clarified the characteristic features of the BgBs interactions [43]. The results appear 
to be reasonable, if we consider i) the spin-fiavor SUq symmetry of Bs, ii) the weak pion effect in the 
strangeness sector, and iii) the effect of the flavor-symmetry breaking (FSB) of the underlying quark 
Hamiltonian. In particular, the existence of the Pauli-forbidden state in the SU3 (ll)s channel (for 
the flavor-symmetric conflguration) and the almost Pauli-forbidden state in the (30) channel (for the 
flavor-antisymmetric conflguration) is of prime importance, in order to understand the origin of the 
Pauli repulsion and the relative strength among various baryon-baryon interactions. This is the most 
prominent feature of the spin-flavor SUq symmetry. The second feature, the weak pion effect in the 
strangeness sector, is a direct consequence of the SU3 relationship in the model framework. The last 
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item is important even in the kinetic-energy term, since the heavier baryon masses reduce the kinetic 
energies appreciably. In fact, the hght pion mass itself is a reflection of the chiral symmetry, which 
gives the most important source of the FSB. 

These B^Bs interactions are now used for the detailed study of the few-body systems, as well as 
medium- weighted hypernuclei and baryonic matter through the G-matrix calculations. In particular, the 
G-matrix calculations are very useful to clarify the characteristics of the model interactions. Both fss2 
and FSS are found to produce the nucleon single-particle (s.p.) potentials and nuclear-matter saturation 
curves, which are similar to those obtained with the standard meson-exchange potentials. For example, 
fss2 gives results very similar to the Bonn model-B potential. It is interesting to note that the deuteron 
properties calculated with fss2 are close to those of model-C, which has a larger D-state probability 
of the deuteron than model-B. Since fss2 reproduces the NN phase shifts at non-relativistic energies 
quite well, the difference of the off-shell effect between our quark model and the other OBEP models 
does not seem to appear prominently, as far as the nuclear saturation properties are concerned. Some 
interesting features of our quark model appear in the predictions for hyperon properties in nuclear 
medium. The A s.p. potential has a depth of 48 MeV (fss2) and 46 MeV (FSS) if the continuous 
prescription is used for intermediate energy spectra. This value is slightly more attractive than the 
value expected from the experimental data of A- hypernuclei [126]. However, the recent analysis of the 
A s.p. potentials in finite nuclei, using the local-density and Thomas-Fermi approximations, seems to 
suggest that this is an appropriate strength of attraction. The E s.p. potential is repulsive, with the 
strength of about 8 MeV (fss2) and 20 MeV (FSS). The origin of this repulsion is due to the strong 
Pauli effect of the (30) component in the EiV(/ = 3/2) ^5*1 state. This result is consistent with the 
recent analysis [127, 131] of the (ii'~,7r^) and {t:~,K~^) experimental data [128, 129, 130]. Future 
experiments are expected to settle the problem of the E s.p. potential and the isospin dependence of 
the EA^ interaction. One of the characteristic features of fss2 is the LS force generated from the scalar- 
meson EMEP. If this contribution is large, the cancellation of the LS and LS^~^ components from the 
FB interaction becomes less prominent. The model fss2 predicts the ratio of the Scheerbaum factors 
Sa/Sn ~ 1/4, which is larger than the FSS value 1/12. Though these ratios arc for symmetric nuclear 
matter with the normal density po = 0.17 fm~^, the density dependence of Sa/Sn is actually weak in 
fss2. 

We also discussed some applications of the A^A^, YN and YY interactions to the Faddccv calculations 
of the three-nucleon bound state, the ANN-TiNN system for the hypcrtriton, the aaA system for ^Bc, 
and the AAa system for aaHg. For these applications, we developed a new three-cluster Faddeev 
formalism which uses two-cluster RGM kernels for the interacting pairs. The model fss2 gives the triton 
binding energy close enough to compare with the experiment. The charge root-mean-square radii of 
and ^He are also correctly reproduced. The application to the hypcrtriton calculation shows that fss2 
gives a result similar to the Nijmegen soft-core model NSC89, except for the appreciable contributions 
of higher partial waves. We find that the ^5*0 state of the AA^ interaction should be only slightly more 
attractive than the ^5*1 state, in order to reproduce the empirical hypertriton binding energy. A desirable 
difference of the ^Sq and ^Si phase shifts around pa — 200 MeV/c is only 0° - 2°, although fss2 is still 
too attractive, giving about 10°. In the application to the aaA system, the aa RGM kernel and some 
appropriate AA^ forces generated from the low-energy phase-shift behavior of fss2 reproduce the ground- 
state and excitation energies of ^Be within 100 - 200 keV accuracy. The very small spin-orbit splitting of 
the excited 5/2+ and 3/2+ states of ^Be is studied in the Faddeev calculation in the jj-coupling scheme, 
where Aa spin-orbit potential is generated from the AA^ Born kernel of the FSS and fss2 RGM kernels. 
It was argued that the P-wavc AA^-SA^ coupling by the LS^"^ and LS^~^a forces is very important to 
achieve the extremely small spin-orbit splitting as an almost complete cancellation between the LS and 
LS^^^ contributions. The EMEP LS force generated from the scalar-meson exchange in fss2 reduces 
this cancellation. As deduced from the NAGARA event [137] for aa^g, the weak AA force is reasonably 
well reproduced by the Faddeev calculation of the AAa system, using the present Aa T-matrix and the 
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full CC AA-SAT-SE T-matrix of fss2. 

The interactions of the octet-baryons, derived from the quark model fss2 in the RGM formalism, 
enable us to reproduce a large number of NN and YN scattering observables. To further test these 
interactions, it is vital to apply them to the systems of more than two particles. The three-cluster 
formalism presented here will serve as a basic tool for this purpose and open a way to solve the few- 
body systems which interact via the quark-model baryon-baryon interactions, while keeping the essential 
features of the RGM kernels, i.e., the non- locality, the energy dependence and a possible existence of the 
pairwise Pauli-forbidden state. Since this formalism makes it possible to relate the underlying NN and 
YN interactions directly to the structure of the few-baryon systems, we expect to learn the interplay 
of the ingredients of the quark model, the baryon-baryon interactions, and the structure and binding 
mechanism of the nuclear systems. 
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Appendix 



Table 18: List of acronyms and abbreviations used in the text 



Term 


Meaning 


Term 


Meaning 




octet-baryons 


QM 


quark model 


NN 


nucleon-nucleon 


FB 


Fermi-Breit 


YN 


hyperon-nucleon 


FSB 


fiavor-symmetry breaking 


YY 


hypcron-hyperon 


CIB 


charge-independence breaking 


BsBg 


baryon-octet baryon-octet 


CSB 


charge-symmetry breaking 


RGM 


resonating-group method 


LS-RGM 


Lippmann-Schwinger RGM 


GCM 


generator-coordinate method 


CC 


coupled-channels 


OCM 


orthogonality condition model 


T-matrix 


RGM T-matrix in Eq. (2.58) 


OPEP 


one-pion exchange potential 


LS 


spin-orbit 


OBEP 


one-boson exchange potential 




antisymmetric LS 


EMEP 


effective-meson exchange potential 




antisymmetric LS with Pg. in Eq. (2.32) 


s.p. 


single-particle 


QLS 


quadratic LS in Eq. (2.33) 


h.o. 


harmonic-oscillator 


is 


one-body LS 
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Table 18: -continued 



Term 


Meaning 


Reference 


KLriVl-r 


quark model including effective-meson exchange 






potentials at the baryon level 


for Qr? o^l 

pO, C)D, 6i\ 


TPQQ 


quark model including scalar- and pseudoscalar- 






meson interactions 


Ton /in /111 

[oy, 4(J, 4iJ 


RGM-H 


an alternative version of FSS 


[40, 41] 


issz 


quark- model including scalar-, pseudoscalar- and 






vector-meson interactions 


[42, 43J 




quark-model package homepage 

http: / / qmpack.homelinux.com/ ~qmpack/index.php 




oAlJJ 


Scattering Analysis Interactive Dial-up 




or yy 


phase shift analysis in the SAID program 


\a^ 1 
[blj 


r WAyo 


phase shift analysis by the Nijmegen group 






Nijmegen soft-core potential 


FlO/l i/lQl 
[iz4, ioo, i4oJ 




Nijmegen extended soft-core potential 


roc o(^l 

[ZO, ZDJ 


*^u-jjonn 


charge-dependent Bonn 




RSC 


Reid soft-core potential 


[160] 


AV18 


Argonne potential with 18 operators including GIB and CSB 


[161] 


SB 


Sparenberg and Baye potential in Eqs. (3.10), (3.12) 


[135, 178] 
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